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DIGITAL COMPUTER ANALYSIS OF CONDENSATION 
IN HIGHLY EXPANDED FLOWS 


SUMMARY 


The application of the IBM 7090 Digital Computer to the theoretical 
prediction of condensation in highly expanded flows is presented. The 
equations of the spontaneous nucleation theory of Frenkel (13) are com- 
bined with the steady one-dimensional diabatic flow equations for a solu- 
tion of the expansion of a pure vapor. A digital computer program for 


the solution of these equations is compiled and presented. 


The theoretical prediction of the condensation of nitrogen is com- 
pared to experimental results and variations in specific heat, latent 
heat of vaporization, surface tension, and rate of expansion are inves- 
tigated. The theoretical calculation is applied to metal vapors and 


the results for copper and zinc vapors are presented. 


The results indicate that the theoretical solution gives a reason--- 
able prediction of the condensation in highly expanded flows. The 
degree of supersaturation increases with an increased rate of expan- 
sion and for a proper set of initial conditions ''condensation free" flow 
is obtainable. The rate of expansion and the surface tension are the 
most critical parameters in the equations for condensing flow. Vari- 
ations of specific heat and latent heat of vaporization show only minor 


effects on the end result. 
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I, INTRODUCTION AND SURVEY OF THE LITERATURE 


The phenomenon of the spontanecus condensation of a gas or vapor 
into the liavid phase has been the subject of numercus theoretical and 
experimental investigations for the past three-auarters of a century. 
The condensation of an expanding vapor was first noted in 1887 by 
R. von Helmholtz; however, this phenomenon did not enter the field of 
aeronautics until about 1940 when its occurrence was noted in the 


supersonic wind tunnel. 


In a supersonic nozzle, the gas or vapor undergoes an isentropic 
expansion which results accordingly in a decrease in the temperature 
and pressure. Generally the saturation vapor pressure decreases with 
temperature more rapidiy than does the static pressure of the isen- 
tropic expansion, hence saturation conditions are approached. If the 
rate of expansicn is very rapid the flow will pass through the point at 
which saturation temperature and pressure are reached without any 
condensate being formed. The vapor is then in a non-equilibrium state, 
supersaturated, in which the existing vapor pressure is higher than the 


corresponding equilibrium saturation pressure. 


The first person to note that a supersaturated vapor can exist was 
R. von Helmholtz, who showed that saturated steam expanding through 
an orifice into the atmosphere would remain clear for some distance 
before becoming cloudy. Stodola (1) furthered this by devising experi- 
ments on superheated steam which would detect the onset of condensa- 


tion and provide eviderce of super saturation. 


Si. 
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In the late 1940's and the early 1950's an intensive series of inves- 
tigations were conducted to determine the condensation of air, water 
vapor in air, and nitrogen. Many of these were conducted at the Cal- 
ifornia Institute of Technology, and worthy of note are the investiga- 
tions of Head (2) on water vapor in air, Willmarth and Nagamatsu (3) 
on nitrogen, and a summary report of GALCIT work by Nagamatsu (4). 
À series of investigations of air in high speed wind tunnels conducted 
during this era include the works of Stever and Rathbun (5), Wegener 


and Reed (6), McLellan (7) and Bogdonoff and Lees (8). 


There was a lag of interest in this area in the middle to late 1950's. 
However, with the advent of the nuclear space propulsion system and 
the use of metallic vapors as the working fluid for hypersonic wind 
tunnels, the interest has been renewed. The Office of Naval Research 
sponsored the recent review of existing theories by Courtney (9) in an 
effort to stimulate this interest. Many of these theories are built 
around the spontaneous self-nucleation concept, and they base the 
collapse of the supersaturated state upon the determination of the rate 
of droplet growth to the critical drop size. The complexity of the 
nucleation phenomenon requires that many simplifying assumptions 
be made in order to evaluate the problem by hand calculations. 
Courtney (9) attempts to obtain a more complete evaluation of the 
nucleation theories by carrying out the calculations on a digital 


computer. 


A recent study of air in hypersonic wind tunnels by F. L. Daum (10) 
indicates that, as yet, there is no proven available theory for predict- 
ing the amount of supersaturation expected in a hypersonic nozzle from 


which may be determined the point of onset of condensation. 
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Until recently none of the investigations had included the problem 
of metal vapors. The investigation by Hill, et al, (11) in 1962 adapted 
the existing theory to the case of metal vapors and provided results 
which indicate that this is a valid extension of the theory. The devel- 
opment of the nucleation equations and flow equations did not prohibit 
their application to metals, but since the need had not existed the 
application was not attempted. The authors noted have not mentioned 
any adaptation of the digital computer to the hypersonic flow problem. 
Hill suggests that this is the next logical step and that a computer pro- 
gram should be used to investigate the effect of variations in surface 
tension, nozzle length, area distribution, inlet conditions, and other 


parameters as well. 


Purpose of the Present Investigation 


At the inception of this study there was no known application of 
condensation theory to the case of metal vapors. The Hypersonic 
Wind Tunnel Group of The University of Michigan was undertaking 
the experimental investigation of the condensation of metal vapors in 
highly expanded flows, and thus arose the need of a method for the 


theoretical prediction of the onset and magnitude of condensing flow. 


Of greater importance was the need for an understanding of the 
effect a variation in any of the various parameters might have on the 
end result. Thereby arose the need for a digital computer program 
from which numerous results can be obtained for multiple combina- 


tions of the various parameters. 


The extent of the present study is to: 
(a) review the literature to determine which of the more promis- 
ing theories will most readily adapt to the IBM 7090 digital 


computer, 





ЕДЕ 


(b) write the program and test against some known experimental 
results to validate the method, 

(c) obtain theoretical results for several metal vapors under 
reasonable hypersonic test conditions, 

(d) vary surface tension, specific heat, heat of vaporization, and 
nozzle parameters to show the magnitude of as on the 


problem. 


The first approach was one toward a general problem with all 
quantities variable and with multicomponent working fluids. This 
rapidly exceeded the capabilities of the computer, and of the author, 
and dictated the simplified approach which is employed. The system 
will be developed for one-dimensional expansion of pure vapors which 
obey the perfect gas law and which have a constant ratio of specific 
heats. Initially the specific heat, latent heat, and surface tension will 


be taken as constant and later will be varied to observe the effect. 





Il THERMODYNAMIC FUNDAMENTALS 


А. Isentropic Expansion of a Vapor 


The assumption that the vapor follows the perfect gas law is, in 
general, well justified and the pressure and temperature can be related 
by 


) Ж Р ў 
£ z T (2.1) 


Any vapor whose latent heat is large, i.e., L > Eu will approach 
saturation as the temperature is decreased by an isentropic expansion, 
and if the expansion is rapid enough there will be a thermal lag so that 


a supersaturated region will be reached. 


The Clausius-Clapeyron equation can be used to relate the variation 
of the saturation vapor pressure with temperature to the latent heat of 
vaporization. 

аға L(T) (3:3 

QT mv vv P 
where p. is the saturation vapor pressure corresponding to a partic- 
ular temperature and for an infinite droplet radius, L(T) is the latent 


heat of vaporization, and ue and v, are the specific volumes of the 


L 
vapor phase and liquid phase, respectively. 


Since the liquid occupies a very small vortion of the total volume, 
| š 1 i 
we will assume Voc аР and Vz D . Assuming L - constant, 


equation (2, 2) becomes И 
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where, = the molecular weight of the vapor 


R = universal gas constant. 


In order to compare the slope of the isentrope with (2. 3) we will 


differentiate (2.1) logarithmically, 


p Sa SO T (2. 4) 
from which, 
С š 
dp "p" p (2. 5) 
dT R T l 


by virtue of y = С DAC and R/u = C - C, . The ratio of slopes is 
therefore 

Po d (2. 6) 

ОП ТС 

р 

Thus for the case of L/T > E the saturation vapor pressure decreases 
with decreasing temperature more rapidly than does the vapor pressure 
in an isentropic expansion, and the supersaturation ratio, p/ Pio? “ontin- 
ually increases until saturation and, eventually, supersaturation are 


obtained. 


Head (2) gives slope ratios for several vapors and shows a typical 
value of 2 for nitrogen nozzles operating at normal reservoir conditions. 
Similarly, he shows that if a vapor had a small latent heat of vaporiza- 
tion, L/T <C_, an isentropic compression is necessary for the approach 


2 


to saturation conditions. 





Upon reaching the saturation point the fluid has a choice of two 
routes by which to continue the expansion. One, the saturated equilib- 
rium expansion, involves a gradual process of condensation which com- 
mences at the saturation point and maintains thermodynamic equilibrium 
throughout the expansion. The other, a supersaturated expansion, is 
characterized by a delay in the onset of condensation to some point down- 
stream of the saturation point at which the vapor temperature is well 


below the corresponding saturation temperature. 


Before continuing further with these alternatives it will be neces- 
sary to develop the diabatic flow equations which account for the effect 


of the condensate. 


B. One-Dimensional Steady Flow Equations With Condensation 


The equations have been derived by many authors; however, we are 
interested here ina specific case which includes the following assump- 


tions: 


1. The vapor phase is a pure vapor and may be treated as a per- 
fect gas. 

2. The condensed mass is a liquid, is uniformly distributed 
throughout the gaseous components, and has the same speed 
and temperature as the stream. 

3. The volume of the condensed phase is negligible compared to 
the total volume. 

4. The nozzle flow is frictionless and 1s without heat transfer. 

5 The saturation curve for the vapor may be adequately approxi- 
mated by the Clausius-Clapeyron equation. 

6. The latent heat of vaporization and the specific heat at constant 


pressure remain constant. 
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Assumption 2 requires some discussion, in that the crystallization 
cf the liquid drops or the sublimation directly to the solid state is prob- 
able at the lower temperatures. In the area concerned, the available 
experimental data on the properties of the vapor and the liquid phases 
is extremely limited, and for the crystalline phase it is practically 
non-existent, Highly extrapolated curves are the only sources, in many 
instances, so the liquid assumption will be as good as the available 
data. The zero drag approximation is supported by Stever (12) in cit- 
ing numerous investigations which verify high fractional condensation 
rates. From the equations of Frenkel (12) it is noted that a small crit- 
ical drop radius is necessary for a high condensation rate. In high 
speed flow the drops are in the nozzle for only a short pericd of time, 
so there is only time for a small increase in size due to growth. Thus, 
for high speed expansions with supersaturation, the drop size will be 
small, and the logical assumption is that the condensate and the vapor 


travel at the same velocity. 


Only the major steps in the development of the flow equations for 
a pure vapor are presented in this section. A more complete develop- 


ment is presented in Wegener (14). 


The continuity equation can be written in its usual form, since the 
volume of the liquid is negligible and the vapor and drops travel at the 
Same velocity. With p taken as the density of the mixture and m as the 


total mass which passes a given location in unit time: 


pUA = m = Constant 
(continuity) p? i OE (2. 7) 


Asa result of the assumptions 2 and 3 the equations of momentum and 


State are: 





(momentum) X = - UdU (2.8) 


dp do ат ав 


977 +r U a 3i (2. 9) 


(state) 


From the first law of thermodynamics, with no external heat sources, 


the energy equation is 


ETE Ыы Жо. (2.10) 


The density of the mixture, p, used in the above equations is 
defined by the equation, p = Dewi PL; The density Pr is the mass 
of liquid condensate per unit volume of the vapor. This requires an 
assumption that the liquid droplets be uniformly distributed throughout 
the vapor so that the densities of the vapor and of the liquid may be 
referenced to the same volume. g is the mass fraction of the working 


fluid which has condensed into the liquid phase. 


By use of the assumption that the vapor and the liquid are at the 


same temperature, and for constant ED and L equation (2. 10) becomes, 


d (U7/2 « C,T - gL) - 0 (2. 11) 


and 


(energy) Бао с ат - L dg = 0 (2112) 


By introducing the sound speed by the perfect gas formula in 


terms of vapor density, A? - y _ T, an alternate form of (2.12) is 


СОС ошаса (2.13) 


(y - 1) Mp + CT 
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The four equations, (2. 7), (2. 8), (2. 9), and (2. 12) involve the six 
unknowns, p, p, T, U, g, and A. For this particular problem, it will 
be assumed that g and A are determined from other sources and p, P, 


T and U are obtained from the above equations. 


Area will be known for a given nozzle, or the area distribution 
must be determined from a given Mach number or pressure distribu- 
tion. The question of determining the condensation, g, depends on the 
assumed flow conditions. If it is assumed that thermodynamic equilib- 
rium exists at all times the Clausius-Clapeyron equation can be used. 
This method will be developed in Section V. However, in the event 
supersaturation takes place, the determination of g is not so simple 


and nucleation theory must provide the missing equation. 


After devoting one section to the determination of the onset of 
condensation, and one to nucleation theory, the flow equations will be 


married to the condensation equations for a complete solution. 





Il, PREDICTION OF THE ONSET OF CONDENSATION 


A review of the literature shows that indeed there is no proven 
available theory which is valid for predicting the onset of condensation 
in a supersaturated flow. A solution of the equations for predicting 
this point always requires an assumption as to the amount of condensate 
formed, the droplet formation rate, or the deviation from isentropic 


flow. 


A marriage of the flow equations (thermodynamic) and the conden- 
sation equation from nucleation theory provide a means of computing 
the flow properties for a given area distribution. A plot of a flow param- 
eter during the period of condensation shows a sharp deviation from a 
plot of the same parameter undergoing isentropic expansion. Experi- 
mental investigations have determined that static pressure, p, and pitot 
pressure, p'o are strongly affected by condensation and a jump in these 


quantities usually is used to detect the onset of condensation. 


Here we are defining the point of onset of condensation as that point 
along the flow at which the initial condensate formed is sufficient to 
cause the computed value of p or T to differ a detectable amount from 


its corresponding value for the isentropic expansion. 


The complete solution could be solved for very small steps starting 
at the saturation point and continuing throughout the length of the nozzle. 
From a plot of actual pressure, p, compared to the pressure distribu- 
tion for an isentropic expansion the point of onset of condensation cou!d 
be determined. However, this would require an excessive number cf 
calculations, and for certain combinations of nozzle parameters and 
stagnation conditions the vapor may traverse the nozzle without ever 


condensing and several minutes of computer time would be wasted. 


Eq 
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Therefore, it would seem expedient to establish a criterion for estimat- 
ing the point of onset of condensation and then to continue the stepwise 


computation of the full set of equations from this point onward. 


Since the expansion of the vapor, in which there is no condensed 
liquid, will be assumed isentropic, any deviation from isentropic con- 
ditions should serve as an indication of condensation. As previously 
stated, the static pressure is sensitive to condensation and this will 
be used as an indicator. From the calculations of Hill (11) the static 
temperature, T, also appears to give a strong indication. At certain 
times it is quite possible that the temperature might give a better nio 


cation than the pressure, and this factor now will be investigated. 


The flow equations from thermodynamics, (2.7), (2.8), (2.9), and 
(2.13), relate isentropic conditions when there is no condensate present, 
g = До = 0, and relate actual conditions when the vapor is condensing. 

A proper combination of these equations should produce a parameter . 
which indicates the deviation from the isentrope. Since we are inter- 


ested in temperature and pressure we will recombine the flow equations 


to obtain expressions for E and > š 
dp. dU dA. 
5 Vis С 0 | (2.7) 
dp = - pUdU (2.8) 
dp do dT dg (2. 9) 


PAI A е 


2dU dT Ldg | 


(y - 1) M UT CT (2.13) 
| p 


Substitute equations (2. 7) and (2. 13) into (2. 9) to obtain 
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Insert this into equation (4. 1) and rearrange 
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| | f 
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Bm МУ. мыш ат 
E 2 FENA OCT e 
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For this particular problem we are interested only in the point of 
onset of condensation. Therefore, we aSsume isentropic conditions 
exist up to the point and the only condensate present will be the dg 
formed over the incremental step at which we are considering the 
onset to take place. The calculations will be accomplished in a step by 
step manner on a digital computer, so we write all increments as A 


and set g = 0 to get 


DES 
Ap __YM ШАН Сн EN (3. 3) 





p M2 - 1 A oe 


j 


From (3. 3) it can be seen that over steps in which there is no conden- 
sate (Ag = 0) the stepwise change in pressure will be isentropic. How- 
ever if Ag £ 0 then the step change in p will differ from isentropic. 


We therefore rewrite (3. 3) as 





=] = 


2 
Ар УМ” АА А І, | 
= x is (3. 4) 
p м2 _ 1 А л (ст 
We now can define 
A L 


as the fractional deviation from isentropic of the pressure change over 
an increment Ax. The magnitude of ер which can be considered a sig- 


nificant deviation is yet to be determined. 


Now we examine a similar analysis of the temperature as a possible 


indicator. From equations (2. 13) and (2. 8) we obtain 


dT Yy-14(.,,9p, U 
те (1 g) + Ç T dg (3. 6) 


` p 





By substituting the value of Ap/p from (3. 4) the equation becomes 








2 
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ee 
T Y M -1 A AA p 
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Tm 2 A “ААС ЕЕ I CE 
ШЕСІ р (y-1)M p 
АТ (у-1) м? АА А Ш ум 
ERE ice e qc Em - Ag (3. 7) 
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Again we define 
TA y-1/M^ p - 1lA (3. 8) 
“Т АА y-1 J Т > 





as the fractional deviation from isentropic of the temperature change 


over AX. 
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Note: A comparison of (3.5)and (3.8) shows a difference only of 


the coefficient of L/C T. Thus as Б 


Er approaches the a of 5 However, an examination of ет апа er 


for given conditions of L/C T, AA, Ag produces: 


ы i.e., as M-—l, 


(a) for condensation in the supersonic portion (M > 1) —— Er P Ж 


(b) for condensation in the subsonic portion (M < 1) — e° > Ep 


As an example we — reasonable values of L/C T z 3and M - 10: 


A 
АА Ag. Therefore, a 
dual criteria will be M. So that ет will determine the onset of 


condensation when M>1 and 5; will be used when M < 1. 


the result gives E - = (2) & x Ag and Єт == (9, 5) — 





IV. NUCLEATION THEORY 


A. Condensation Nuclei 


In the process of condensation, the saturated state is referenced 
to a plane liquid phase and the vapor is considered to condense ona 
flat liquid surface or ona "'cold" container wall such as the moisture 
formed on the outside surface of a glass of ice water. However, in 
highly expanded flows the liquid surface is not present and the "cold" 
wall is not available, due either to the "hot" boundary layer surrounding 
the flow in a nozzle or to the total absence of a wall in the case of the 


free jet. 


The vapor then must look for other surfaces on which to condense. 
These surfaces appear either in the form of impurities, such as par- 
ticles of dust, etc. , or as small drops formed of clusters of molecules 


of the vapor itself united by statistical collision. 


Obviously, for a truly pure vapor there would be no foreign nuclei; 
however, dust particles in clean air have been estimated upward of 103 
particles per cubic centimeter, so this possibility must be considered. 
Head (2) has shown that these foreign nuclei do play an important role 
in the onset of condensation. However, for highly expanded flows 
where the stay time is very short, Stever (12) has shown that the high 
rate of condensation could not be caused by condensation onto the for- 
eign nuclei alone. Stever (12) quotes an example from Oswatitsch's 
study of water vapor which considers 10? foreign nuclei per cubic cen- 
timeter traveling in a gas for a distance of 10 centimeters where there 
is super cooling of 30°С. А slightly supersonic velocity of 3.3 x 10° 
centimeters per second would traverse this 10 centimeters in 3 x 1074 


seconds. Using the faster growth formula for droplet growth gives a 
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radius of growth of about 3 x joie centimeters in this distance. The 
amount of water condensed over this span then is 1072 grams per cubic 
centimeter, which releases an amount of heat that is very slight and 
which has an insignificant effect on the flow quantities. The actual 
condensation process has a marked effect upon the flow properties; 
therefore it must be concluded that the number of drops far exceeds 
the number of foreign nuclei. We must then resort to the spontaneous 


formation of the droplets onto which condensate growth occurs. 


There is an opening in the above argument, in that the growth rate 
onto the foreign nuclei could be grossly in error. However, Head (2) 
found droplet formation rates on the order of 101° drops per cubic 
centimeter per second in his investigation of water vapor; thus, the 


high particle formation rate tends to support Stever's argument. 


For this analysis it will be assumed that the condensation phenome- 
non is one of spontaneously formed nuclei and is triggered by the attain- 
ment of critical size droplets onto which the supersaturated vapor 


condenses. 


The concept of the critical drop size as developed in Stever (12) 
stems from the energy required to evaporate a drop of liquid which 


contains n molecules. 
2 
E(n) =nE, - 4nró0 (4, 1) 


In (4,1) E, is the energy required to remove by evaporation a single 
molecule from an infinite surface of liquid and o is the surface tension. 


To evaporate one molecule would require 


= _ liq 
AE(n) = E y 3x (4, 2) 
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where V is the volume of a molecule of the liquid. An application 


liq 
of Boltzmann's law gives vapor pressure 
_ АЕ 
оК KT 


and the equilibrium vapor pressure for a plane surface 


Combining these gives the Thompson-Gibbs equation 


2У.. o/rkT 
De і (4, 3) 


Poo 


Expressing the molecular volume Ya as WN Dy and rewriting it as 


a log function, (4. 3) becomes 


aug 20 


m р/р, x rp, kN,T y rp, RT 





(4. 4) 


Note that with the occurrence of r in the denominator of equation 
(4. 4) the equilibrium between the vapor and the liquid drop is unstable. 
For a vapor which is saturated with respect to a given drcp size, if 
one molecule sticks to the drop the radius increases and the existing 
vapor pressure is greater than the new equilibrium vapor pressure; 
thus, condensation will continue. The reverse is true for evaporation, 
in that if one molecule evaporates the existing pressure is lower than 


equilibrium and evaporation continues. 


Thus arises the concept of the critical drop radius, r*, for which 
r < r* the droplet vaporizes and for r » r* the droplet grows. Solving 


(4, 4) gives 
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r* = 

Hill (11) points out that in a constant environment, if the critical 
radius was used to specify the initial drop size then the growth rate 
would be zero and there would be no condensation because at critical 
size the decay and growth are equally probable. Hill bases his use 
of 1. 3 times the critical radius on Becker and Doering's conclusion 
that at that size the probability of decay is close to zero. In this anal- 
ysis the vapor will be flowing and any critical drops formed will not 
remain in a constant environment very long since the quasi-steady 
steps for the calculations will be very small. Therefore, the com- 


monly used r* of equation (4, 5) will be used as the initial drop size. 


B. Kinetics of Nucleation and Condensation 


As implied previously, the spontaneous nucleation approach will 
be the only one considered in this analysis. There have been volumes 
written in the development of this area and no attempt at reproduction 
will be made here. The general approach and final result of the most 


frequently referenced authors will be stated, however. 


In highly expanded flows the state of the vapor can change very 
rapidly from the unsaturated state to the highly supersaturated state. 
Frenkel (13) treats the equilibrium distribution of drops in an unsatu- 
rated state as a dilute solution of different substances in the vapor as 
a solvent. By classifying the different solutes by the number of mole- 
cules, n, they have in a drop, he obtained the expression for the 
number of drops containing n molecules, 

1 2 
Ж a = 
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where N - number of molecules in the vapor state 


‚Ф 


ES фр” = thermodynamic potential of single molecule 


4т E - surface energy term. 


This droplet distribution is based on the assumption of equilibrium 
< Ф 


tion can be used for supersaturated conditions (4 


); however, it is generally assumed that the samé distribu- 
a = $): Note that 
for the supersaturated condition the number of drops increases with 
increasing radius and the distribution is unstable. Thus, it becomes 
essential to determine the rate at which the distribution changes from 
the stable distribution of the unsaturated state to the unstable distribu- 


tion of the supersaturated state. 


This requires a determination of the rate of formation of critical 


sized drops, J, in drops per cubic centimeter per second. 


Head (2) states the equilibrium equation of Volmer 


AW 4T r* 2 


E ^. 8kT 


J=K-e K:e (4, 7) 


where AW is the total work of the formation of the droplet. 


The determination of the constant, K, has been developed by sev- 
eral authors, including Volmer, Becker and Doring, and Frenkel, 
All assume a quasi-stationary phase transition process in which the 
number of molecules which condense are instantly replaced in the 
vapor such that the number of vapor molecules is maintained constant. 
Volmer considered only the number of drops reaching critica! size, 
the number of vapor molecules striking the surface of the critical size 
drop per unit time, and the number of molecules in a drop of critical 


size. 
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Becker and Doring derive the nuclei formation frequency purely 
from kinetic considerations and take into account the evaporation from 
the surface of the drop, as well as the condensation onto it. This ap- 
proach is considered by most authors to be better than others; however 


their solution is rather difficult to apply. 


Frenkel uses the method of Zeldovich to integrate the differential 
equations of Becker and Doring and gets a much simpler derivation of 
an equation that is almost identical to that of Becker and Doring. The 


ratio of the result of Frenkel to that of Becker-Doring is (п) : and 
this is considered minor. n* is the number of molecules in a critical 


sized drop. 


A partial derivation of the nucleation rate equation is presented 
in Appendix A. This derivation is necessitated by an apparent error 


of 7 in the original equation presented by Frenkel (13). 


For this analysis, equation (A. 7) will be applied as the nucleation 
rate equation. 
2 _ 4r ers 
| p | zit 210 к 3kT 
P LS 





(4. 8) 


where J = number of critical drops formed per unit volume per unit 


time. 


C. Limitations 


1. Generally o as a function of temperature is not known and this 
is a large factor in the result. Surface tension appears to 
the third power in the exponent; hence a small error will have 


a large effect upon J. 
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2. The rapidity of expansion will show marked deviations from 
the assumed quasi-stationary condition. 
3. The condensing flow may have sublimation directly from vapor 
to solid which is not considered. 
4. At high degrees of supersaturation the drop radius will be small, 
and Head (2) has shown that the variation of o with r is large 
for drops of few molecules. Stever and Rathbun (5) have included 


this effect into the Frenkel equation to obtain 
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Bogdonoff and Lees (8) have an approach which is entirely dif- 
ferent from that of Head and Stever and which produces quite 
different results. Until these theories have been developed 
further, it is not practical to include either theorem in this 


analysis. 


From the previous discussion of the assumptions made in the formu- 
lation of the nucleation theory and of the limitations above, it is apparent 
that the best one can expect in an analysis such as this will be an order 
of magnitude estimate of the condensation parameters and a qualitative 


analysis of the flow conditions. 


D. Growth Estimates 


Heretofore, in order to calculate the formation of the critical drops, 
we have considered the temperature of the liquid drops to be the same 


as the surrounding vapor. In actuality, the impinging and vapcrizing of 
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molecules onto and away from the droplet leave a surface temperature 
that approaches the saturated equilibrium value for the surrounding 
vapor. This temperature difference allows for growth of the particle. 
The growth of small drops, i.e. , r «« mean free path, occurs by the 
molecular transport of molecules onto the surface, some of which con- 
dense, and by the conduction away of latent heat of vaporization by the 
evaporation of a fewer number of molecules which leave at higher kinetic 


energies. 


Oswatitsch laid the basic groundwork for growth in a flowing vapor. 
Stever (12) follows the method of Oswatitsch to develop an equation for 


the growth rate of the small drop. 


The development presented by Stever makes use of the kinetic 
theory of free molecule flow (See Patterson (15)). Consider a drop 


of radius r, temperature T,. , and saturation pressure (ps, ) . with 


Па Па 
a surrounding vapor of pressure p and temperature T. The heat trans- 


fer from drop to vapor per unit time is 


k 
ais er a US 


where o is an accommodation coefficient which varies from 0 to 1, 


depending on the vapor and liguid surface. 


Also from kinetic theory the heat transferred to the drop surface 
due to condensation is Lfp /KT/2r m where p VET/27 m is the amount 
of vapor (in the form of molecules) impinging on a unit surface and f is 
the fraction which condenses. The energy being transferred to the sur- 
face is equated to the energy being transported away to obtain an expres- 


sion for f; 
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к 2a (Tig - T) 
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the mass Mp of the drop is then increased by 


uie a Er a eee 


S Py and dM, = 47 s Рі. dr, 


If the calculation is made in quasi-steady steps, time can be expressed 


as a function of step length and 


Аг ya LE v ОТ nem Ча) 


where the accommodation coefficient, @, is a constant, approximately 


equal to 1 and T,. is the temperature of the surface of the liquid drop. 


liq 
The accommodation coefficient o is expressed by Kennard (16) as 
а (Е, - Е) = (E, - Е) 


where E. = Energy brought up to a unit area per unit time by the 
impinging molecules 
ES - Energy carried away from the molecules leaving a 
unit surface per unit time 
E - Energy that the departing molecules would carry away 
if they carried away the same mean energy per mole- 
cule as does a stream issuing from a gas in equilibrium 


at the surface temperature. 
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Patterson (15) lists coefficients of the range from .8 to 1.0, 
and Stalder (17) found that for free molecule flow an accommodation 
coefficient of 0.9 checks well with test results. Thus we arrive at 
the statement, "approximately equal to 1. " This only applies for cases 


where free molecule flow can be assumed to prevail. 


On the other hand, where r approaches the magnitude of the mean 
free path the heat removal is controlled by normal gaseous heat con- 
duction. For this large drop Stever (12) gives the growth equation for 
period of time, t, 

z (Tuo ^ D 
NE 


Other authors have made similar analyses, and Wegener (14) re- 


Ar = і. (4. 11) 


cords the quasi-steady solution of Buhler (18) which includes an expres- 


sion for the condensation coefficient. 


1 
2 (T. - T) 


20 
! y rp, b 


where T the surface temperature of the drop, may be chosen as the 


1 


ny (4, 12) 











equilibrium saturation temperature of the vapor with respect to the 
drop. As in the previous case, Wegener states (4. 12) is valid only 
for small drops and the equation recommended for large drops is 


merely a revised form of (4. 11). 


Hill (11) notes that typical calculations of condensation nuclei indi- 
cate very small sizes which are generally much smaller than the mean 


free paths; thus only the free molecule equations need be considered. 
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His analysis produces a growth rate expression of 


dr _ 3p [RT 


2 кале - RT] (4. 13) 
en 20, PL VRT 9 


in which Ve = change of internal energy due to condensation. The 
term Ту is droplet temperature, ard Hil! presents calculations to show 
that the rate of drop temperature change is enormous compared to the 
rate of growth. For water vapor the temperature approaches wet bulb 
temperature in approximately the same time that it takes the number 
of mclecuies in the drop to increase by a factor of 10. Hill uses this 

in expanation of the assumption used by Oswatitsch that the drops al- 


ways have wet bulb temperature during the growth period. 


After a displacement, the drop temperature will rapidly approach 
the temperature at which d T dt = 0 at a much greater rate of change 


than dr/dt. Hill uses this to form an expression for T_ which depends 


D 
only upon the instantaneous values of r, p and T. 


2 Yig 1-29 exp 22.29 е (4 14) 
3 RT р p, RTyr Tp T | 


where & is the fraction of impinging molecules which condense. 





From the investigations on water vapor, nitrogen, air and a few 
metal vapors it appears that for highly expanded flows the small drop 
size is the general rule. The calculations will then be initiated by 
applying (4.10) as the growth rate equation and choosing drop tempera- 
ture as the equilibrium saturation temperature of the vapor. The ac- 
commodation coefficient in (4. 10) will be varied to study the effect and 
a will be placed into the computer program as a data input so that accu- 


rate values, when known, can be used for the calculation. 





V. ANALYSIS 


The basic equations are established in Sections H through IV. The 
proper combination of these equations produces a theoretical solution 
for the hypersonic expansion of a pure vapor in which a supersaturated 
state occurs. Another solution is obtained for which there is saturated 
equilibrium expansion in lieu of the supersaturated portion. Both solu- 


tions, however, stem from the same basic model. 


The basic model considered here consists of assuming: 

1. ‘One-dimensional nozzle" with geometrically specified inlet 
and diffuser. 

2. A pure vapor with specified stagnation conditions and known: 
constant specific heats, latent heat of vaporization; surface 
tension, for infinite radius as a function of temperature; and 
saturation curve. 

3. An isentropic expansion to the point of onset of condensation 
and then diabatic flow through the region in which condensa- 


tion is taking place. 


A. Basic Assumptions 
In addition to the assumptions of the basic model, this analysis 
also embodies the multitude of simplifying assumptions which have 
been included in the derivation of the basic equations. The type of 
equations and the specific assumptions made are: 
1. Flow equations 
Mass flow is constant 
Flow is one-dimensional and steady state 
c. Volume of condensed phase is negligible when compared 


to the total volume 
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e. 


28 


Nozzle is frictionless with no heat transfer across nozzle 
wall. 

Condensed mass is liquid, is uniformly distributed through- 
out the gaseous components, and has same speed and tem- 


perature as the stream 


Prediction of onset of condensation equations 


a. 
b. 


This embodies same assumptions as flow equations 
No condensate is formed prior to the predicted point of 


onset of condensation 


Nucleation equations 


а. 
b. 


d. 


е, 


Drops are formed by spontaneous nucleation only 

Orly drops reaching critical drop size continue to grow 
Number of molecules of vapor is maintained constant 
Saturation pressure corresponds to saturation pressure 
of a droplet of infinite radius 

The ordered velocity of the particles is neglected 
Equilibrium values of the probabilities that one molecule 
will leave from or cordense on a unit surface and equilib- 
rium pariicle distributions are assumed valid for the 


non-equilibrium case 


Growth rate equations 


Drop radius is much less than mean free path 
Accommodation coefficient, o, is approximately equal to 
1.0 

Drop temperature is assumed equal to the saturation tem- 


perature that corresponds to the pressure of the vapor 
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B. Supersaturated Expansion 


The supersaturated expansion passes the saturation point without 
any condensate being formed; therefore, a four step solution is required: 
1. Isentropic equations from Shapiro (19) and NACA Report 1135, 


are used to calculate the expansion of the vapor to the point of 























condensation. 
Y 
y-1 
T 
P - P, [| (5. 1) 
a 
1 
| 2 E Y 
М = ЛЫ ПЫ. (5. 2) 
ИКЕ 
y- 1 
e (5.3) 
y+1 
| 2 (y - 1) 
_A*|_2 el 2) 
A-zril* ; M (5. 4) 
1 
2 
U=[2C (T - T)] (5. 5) 


At any particular value of T the above equations are readily 
evaluated for given Po То Dy Y, E and A* The station at 
which these values apply is then determined from the given noz- 
zle geometry: 

2. The saturation point is established by an iterative balance of 


the saturation curve equation with the isentropic flow equation. 
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The saturation curve can generally be expressed in the 
form In р = В - = , in which B and C are constants whose 
value depend upon the choice of units for p and T. The equa- 


tions which require joint solution are 


C 





-— n wm (5. 6) 
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The isentropic relations above are then utilized to calcu- 


late M P T 


sat’ sat’ sat’ A sat P sat? and Un 


E 
The point of onset of condensation is determined by taking 
steps, AT, of supersaturation and computing the flow param- 
eters from the isentropic relations for each successive step. 
The isentropic values of each incremental step in temperature 
are introduced into the nucleation equations for critical drop 


size and formation rate. 
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r* = —— U (4, 5) 
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Note that choice of a very small AT will cause the first 
step to remain very near the saturation point and the resultant 
р/р, = 1 will cause r* to approach an infinite radius. r* falls 
off very rapidly, however, as the degree of supersaturation 


increases, and for most cases a choice of AT > 5k is sufficient. 
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The growth rate is zero for the critical droplet at the instant 


of formation; therefore the nucleation rate for this increment is 


Ni. =: Ir An Ах. (5,8) 


The mass fraction of condensate formed in this increment is 


ы К 3 


"mE * 
T 3m `T T (5. 9) 


Ag 


The epsilon equations from the onset of condensation develop- 
ment, (3. 5) and (3. 8), are now used to determine if this amount 
of condensate is sufficient to cause an appreciable deviation 


from the isentropic expansion. 
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If the € is less than the assumed critical value then the same 
calculations are carried out at step T + AT until the critical 
value of € is equalled or exceeded. If exceeded, a bracketing 


procedure is employed to improve the estimate. 


All quantities computed for Ag at values ofe<e_... 
critical 
are discarded, and it is assumed that there has been no con- 


densate formed prior toes e ... ,. 
critical 


The choice of Ax in equation (5. 8) will effect directly the 
magnitude of N and of Ag in equation (5.9). AA in equations 
(3. 5) and (3. 8) is also a function of Ax so that in the evalua- 


tion of € the effect of Ax is nullified. Calculations have been 
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performed for values of Ax ranging from . 01 to 1.0 cm to verify 
that the quantity has no effect on the prediction of the point of 
onset of condensation. Since the Ax that will be used for incre- 
menting the calculations through the condensing portion of the 
nozzle will be a required input to the computer program, it is 


advantageous to use that same Ax in evaluating (5. 8). 


The condensing portion of the flow requires a joint solution of 
the nucleation and growth equations and the diabatic flow equa- 
tions. The calculations are performed for increments of Ax 
starting from the onset of condensation. Values are estimated 


of p, T and U at the point i and Ag. is calculated from the 
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The subscripts i and j are numbered increments of x 
starting from the condensation point. i acts as a label for 
each group of droplets of a particular size and denotes the 
particular increment in which these drops originated as crit- 
ical sized drops. j denotes the increments presently under 
consideration. Thus a designation r,, denotes droplets 
which were formed as critical drops in increment 1, have 
undergone the growth Ar, and are now undergoing the 
growth Aro. 
Using the above approximate values of Ag and g and the 


unknown values of AA and A the equations from Section II 


are solved for p, T, p and U. 





Ap, AU, AA 
5 der u (2 7) 
Ар | UAU _ (2. 8a) 
p R 
l mh ` T 
( 8) 
Ар Ар АТ Ар; ТЕ 
Ар ар Al 2,9 
p p T ess) Т 
UAU + С, AT - LAg =0 (2. 12) 


The solution of the flow equations is accomplished by the 
following iterative procedure: 

a. assume AU and compute т = ы 217 АТ 

b. solve (2.12) for AT and compute T, = Т. NE AT 

solve (2. 7) for Ap/p 

d. solve (2.9) for Ap/p 

e. solve (2.8a) for AU and check with value from step (a). 
Improve the assumed AU and continue until a solution is 


obtained. 
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These corrected values of p, T and U are now used to 
resolve the equations for Ag, and the process is continued 
until a set of flow parameters which satisfy both sets of 


equations is obtained. 


This iterative procedure is continued for each increment 


Ax up to a specified X nax’ usually the end of the nozzle. 


This set of equations could be evaluated from the satura- 
tion point onward in order to obtain a more accurate predic- 
tion of the onset of condensation. It is obvious, however, that 
the magnitude of the calculations involved is formidable; and, 
the slight error induced by attempting to predict the point of 
onset of condensation, does not justify the refinement. This 
is especially true since the approximations in the basic 


equations are so gross. 


C. Saturated Equilibrium Expansion 


The saturated equilibrium expansion is a gradual process with 
condensation starting at the saturation point and equilibrium flow 
being maintained throughout. This occasionally occurs for slow 
expansions of certain vapors and for cases in which the impurity 
content is high. In this analysis the saturated equilibrium expansion 


is presented only as a comparison with the supersaturated case. 


The calculations are greatly simplified in that the mass fraction 
of condensate can now be specified by the Clausius-Clapeyron 


equation 


dp Бир (2. 3) 





-35- 


This equation and the flow equations (2. 7), (2. 8), (2. 9) and (2. 12) 
usually are integrated and expressions are written for each param- 
eter (see Wegener (14) or Daum (10)). By proper manipulation and 
combination the integrated equations can be reduced to one involving 
g. 
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отт ыт 
sat 
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However, in this analysis the flow equations are programmed in the 
incremental form for the supersaturated case, so it is convenient to 
write (5. 13) in the form 


С T. 


e AS ЕЕ 
sat 


р (5. 14) 








Т | 51-1 
sat 








and to use this for the Ag expression in a solution similar to that 


used for supersaturation. 


1. The saturation point is established in the same manner as 
for the supersaturated case. Small increments of length, 
Ax, are started at the saturated point. 

2. An incremental increase in the mass fraction of condensate, 
Ag, is assumed for the step Ax. With assumed value of Ag 


and known value of A the flow equations are solved for p, p, 


and T. 
Ap AU AA_ 9 
u aA (2.7) 
Ap UAU 


ы (2. 8а) 
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р р Т (1-6) 





(2. 9) 
UAU + C, AT - LAg - 0 (2. 12) 


Equation (5. 14) is solved to check and correct the assumed value of 
Ag and an iterative procedure is established to obtain a satisfactory 
solution of the five equations. Again this is continued in a step by 


step process until the solution is complete.. 





VI. DIGITAL COMPUTER PROGRAM 


A major task in this study has been the adaptation of the mathematical 
analysis of Section V to the digital computer. The description contained 
herein is presented for the future user of the program and in no way 
covers the myriad of detours encountered between the input and output 


statements. 


The Computer Center of The University of Michigan has developed 
a Michigan Algorithmic Decoder (MAD) computer program which com- 
piles much faster than the commonly used FORTRAN program. The 
computer program for this investigation has been written in the MAD 
language and can be translated into FORTRAN if desired. The IBM 7090 
computer program for the condensation of a pure vapor is presented in 


Appendix B. 


The program has been subdivided into a main program and six 
Subroutines. The subroutines enable easier modification of the program 
and provide a greater flexibility in its use. Statements as to the func- 
tion of each of the calculations, the input parameters, and the obtainable 
results are provided in the program at the beginning of each subroutine. 
Appendix B also includes a list of the program symbolic names assigned 


to the parameters used in the equations of Section V. 
A brief summary of the subroutines follows: 


1. Isentropic flows (IFLOWS)—Given T, this program solves 
the isentropic flow equations for p, p, U, M, and A. 

2. Condensation Flow (CFLOW)— For given g and A, this pro- 
gram solves the diabatic flow equations (2. 7 through 2. 12) 
for p, T, p, M, and U. 


ЕТЕ 
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Nucleation Theory I (NUCLE1)— From nucleation theory equa- 
tions (4.5, 4.8, 5.8 and 5.9), NUCLE1 computes the condensate 
formed due to the nucleation of critical sized drops only, i.e., 
growth is not included in this portion. This subroutine is used 
to determine condensate formed for that portion of the main 
program which determines the point of onset of condensation. 
Nucleation Theory II (NUCLE2)— This subroutine evaluates 
the growth equation (4. 10). NUCLE2 is used in conjunction 
with NUCLEI to determine the amount of condensate present 
for all calculations downstream from the point of onset of 
condensation. 

Nozzle (NOZZLE)— The geometry of the nozzle is placed in 
this subroutine, from which A, AA, or x can be determined. 
Vapor Parameters (FRHOL, FL, FSIGMA)—The parameters 
Py: L and o are placed in this subroutine as functions of T 

or as constants and are evaluated as called upon by the 


Main Program. 


The Main Program consists primarily of a proper utilization of 


the subroutines. The operations by the Main Program appear in the 


following order: 


1. 
2. 


Read in the input data. 

Use IFLOWS to expand isentropically from stagnation condi- 
tions. 

Match isentropic equations to saturation curve to determine 
saturation point. 

Use IFLOWS for expansion past the saturation point and 
calculate Ag from NUCLEI. 

Evaluate epsilon equations (3. 5 or 3. 8) to determine onset 


of condensation. 





esp: 


6. Downstream of the onset of condensation, a joint solution of 
NUCLE1, NUCLE2, and CFLOW is obtained at each incre- 
ment, Ax, and the parameters x, N, Ax, A, AA, P, p, T, 


UPM ES, ги, Ar, Ag and g are stored for later print out. 


The saturated equilibrium expansion is compiled in a very simple 


program and need not be presented at this time. 


The Nozzle and Vapor Parameter Subroutines are included ina 
Program Common such that any of the included parameters, when used 


in the Main Program, are automatically evaluated and returned. 


Conditional statements are inserted in the program to account 
for the possibility of condensation occurring in either the converging 
or diverging portions of the nozzle, the possibility of reaching the end 
of the nozzle without any condensation having occurred, and the failure 


of the iteration of NUCLE1 and NUCLE2 with CFLOW to converge. 


The development of the equations is such that any consistent set 
of units can be used for the problem. However, it was found that 
minimum scaling was necessary if the problem was worked in a cgs 
system of units. The problem, as written, is compatible with the 
cgs system of units and a list of the proper input data is presented 


in Appendix C. 





VI. CALCULATIONS AND RESULTS 


The calculations performed in this study are divided into four 


parts: 


1. 


A test of the program by hand calculations and a test of the 
method by comparing the results for a nitrogen expansion 
with the experimental results of Willmarth (3). 

Justification of the assumptions of small radii, no drag, and 
accommodation coefficient approximately equal to one in the 
basic assumptions. Verification of the prediction procedure 
for the onset of condensation. 

A variation of latent heat, specific heat, surface tension, and 
nozzle geometry for the nitrogen expansion to observe the 
induced effect of each parameter. 

Application of the theory to metal vapors by computing the 
expansion from an initial pressure of approximately one 


atmosphere of copper, tin, lead, and zinc vapors. 


A. Validation of the Program 


The input data for the calculations performed in this study are 


tabulated in Table C-1 of Appendix C. A reference to Nitrogen (1), 


Zinc (1), etc. , indicates the set of input data from which the indicated 


results were obtained. Any additional modifications of the input data, 


such as a change ino, L, or nozzle geometry, are listed directly on 


the plot of results. An explanation of the input data required, the 


proper units for each parameter, and the references (20-24) from 


which the vapor properties were obtained is also presented in 


Appendix С. 


-40- 





-41- 


Initially the computer program was divided into small segments 
and hand calculations were performed for Nitrogen (I) to verify the 
computer results. With the assurance that each small segment was 
working properly, the entire program was assembled and hand calcu- 
lations were performed to spot check the calculations of the saturation 
point, the point of onset of condensation, and the first few incremental 


steps of the condensing flow regime. 


With this degree of confidence in íhe program, the next step was 
to check the theoretical results with known experimental results. Тһе 
Nitrogen (1) input data is designed to approximate an experimental 


test of the condensation of nitrogen performed by Willmarth (3). 


Nitrogen properties from NBS Circular 564 (24) and the Interna- 
tional Critical Tables (23) required extrapolation below 70°K for the 
surface tension and the latent heat of vaporization. These values are 
generally valid only to within +10% and, as will be shown later, this 
becomes extremely critical in the case of the surface tension. The 
vapor pressure curve for nitrogen is reasonably well known, and it 


was approximated by the Clausius-Clapeyron equation (see Fig. C-1). 


For comparison with experiment the calculations for nitrogen were 
carried out in full from the saturation point onward and no attempt was 
made to predict the point of onset of condensation. A comparison of 
the experimental results of Willmarth (3) to the theoretical prediction 
of the condensation of Nitrogen (I) is presented in Figure 1. The the- 
oretical curve for Nitrogen (I) shows a degree of supersaturation of 
approximately 3. 5 degrees (20%) lower than is indicated by experimental 
results. This is expected in that the amount of impurities in the experi- 
mental nitrogen was unknown and this advanced nucleation would cause 


an early break. However, of greater importance is the degree of 
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uncertainty in the extrapolation of the available surface tension data. 
In anticipation of the effect of o on the theoretical results, a second 
extrapolation was chosen for o. Theo equation for Nitrogen (II) re- 
mains within the + 10% accuracy of the original surface tension data 
and, in actuality, gives a value for o which is only 7% less than the 


o from Nitrogen (I) at the onset of condensation. 


The theoretical prediction for the condensation of Nitrogen (II) 
compares most favorably with the experimental results in Figure 1. 
The strong influence shown by only a 7% change in o will be discussed 
more fully at a later point in the discussion. Figure 2 presents a 
more complete comparison of the theoretical results for Nitrogen (II) 
as compared to the experimental results of Willmarth (3). The isen- 
tropic and saturated equilibrium expansions are also presented in this 
figure. The theoretically calculated results presented in Figure 2 
employ the full calculation of the nucleation and flow equations for the 
entire length of the nozzle. This is significant in that the flow prop- 
erties prior to the ''condensation shock" lie on the isentrope and the 
flow properties after the ''shock" tend to approach the saturated 
equilibrium expansion. This is the expected result, and it gives 


further confidence in the program. 


The experimental points of Willmarth (3) fall on the isentrope as 
well as along the "condensation shock", which indicates a proper esti- 
mate of the physical expansion and validates the program. The close 
correlation between theoretical and experimental results within the 
"shock" itself lends further support to the basic nucleation theory of 


the condensation process. 





direct measurement. The effective area ratio, (A/A*) 


23. 


Willmarth (3) obtained the pressure results, p, Po and ре! ру 


eff? tempera- 


ture, T, and mass fraction of condensate, g, were then calculated 


from the measured quantities by a joint solution of the flow equations. 


B. 


Justification of Assumptions 


1. Drop Size 

A typical distribution of maximum and minimum droplet 
radius is presented in Figure 3 for the Nitrogen (I) expansion. 
Note the very short span over which the nucleation of critical 
sized drops, J curve, takes place. The average critical sized 
drop has a radius of 4.5 x a cm, and the largest particles 
which leave the nozzle have a radius of 8 x 107! cm. The radius 
of the nitrogen molecule from Loeb (25) is approximately 
Leo" x 107% cm and calculated mean free paths vary from the 
order of We cm at the onset of condensation to Т сіп пеаг 
the nozzle exit. According to Stever (5) the drop radius is 
larger than the minimum required to contain the 10-12 molecules 


necessary for the application of macroscopic theory. 


Therefore it is not unreasonable to assume that small drops 
are the general rule, that the no drag approximation is valid, 
and that since r << mean free path the samll drop growth equa- 


tion can be used. 


2. Accommodation Coefficient 

The effect of a variation in q is presented in Figure 4. Small 
variations near the value of a = 1. 0 have very little effect. How- 
ever, since the growth equation is a direct function of o, the curve 
for a= 0.1 differs markedly from @=1.0. Fortunately, the known 


values for o are from 0. 8 to 1.0 so that the choice of a=1.0isa 
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reasonable one. The computer program is so constructed that 
accurate values of q, when Known, can be inserted as input 
data. 

An additional feature of Figure 4 is the indication that growth 
has little effect up to the "break" point but thereafter controls the 
remainder of the condensation process. This is further brought 
out by the plot of J in Figure 3, which indicates the rapid nuclea- 
tion of a large number of particles to start the "shock'' and then 
a continuing growth of these particles until saturated equilibrium 


conditions are met. 
3. Point of Onset of Condensation 


From the previous discussion and from Figure 3, it appears 
that little or nothing happens in regard to condensation until a 
degree of supersaturation is reached for which the critical drop 
radius is very small. As r* from equation (4. 5) becomes small 
the exponential term, which contains po ceases to be the pre- 


dominant term in equation (4. 8) and J increases rapidly. 


| "ier 
ud (4. 5) 
p, RT In p/p 


9 : 4T ork? 
Та E 2 3-е ЗЕТ (4. 8) 
Pr, A 


The epsilon of Section II should then predict the approximate 
point of onset of condensation. Several calculations for various 
values of € are presented in Figure 5. є = 0 indicates that the full 
calculation of the nucleation and flow equations was carried out 


from the saturation point onward. The choice of epsilon results 
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in only a minor error when compared to the possible 10% error 


I nel“ give almost 


in surface tension. Epsilon values of 10° 
identical results, thus e =. 001 appears to be a suitable choice 
for the expansion of nitrogen. 

Also in agreement are Head's (2) calculations for water vapor 
in which he indicates that the attainment of J = 1)“ droplets per 
centimeter cubed per second signifies the point at which the con- 
densation shocx will break away from the isentrope. Several 
calculations on the computer indicate that J = Tone at the breax- 
away point for nitrogen also. At present, however, this author 
sees no justification for predicting that this will be true for other 


vapors as well. 


Variation of Parameters 


1. Surface Tension 

As indicated previously, the variation of o is the most critical 
aspect of the theoretical prediction of the condensation phenomenon. 
From equations (4. 5) and (4. 8) it can be seen that o enters the 
exponential term of J as a third power. Thus for a 10% change 
in o the nucleation rate changes by several orders of magnitude. 
Figure 6 and Figure 7 present several variations in o to show the 
large effect that this parameter has on the degree of supersatura- 
tion. Later results will show that the effect of other parameters 


is minor when compred to surface tension. 


There appears here a possible method of improving the surface 
tension data. The results of a closely controlled experimental vapor 
expansion could be used aS a reference plot and surface tension 


varied on the computer until a matching plot is cbtained. Providing 
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the other parameters are reasonably well known, this should pro- 
vide a better estimate of o than merely extrapolating a curve far 
away from the few experimentally determined values that are 


available. 


2. Specific Heat | 

сз is reasonably well known and only varies a few percent 
over the range considered in this analysis. For this analysis 
= is chosen as a constant, and from the plot in Figure 8 the 
constant assumption for КЕ appears reasonable. The effect of 
a 10% variation in G is negligible when compared to an equiva- 


lent uncertainty in o. 


3. Latent Heat of Vaporization 

Latent heat is generally a weak function of temperature. For 
Nitrogen (I) the variation in L is only 9% over the range of tem- 
perature from saturation to the nozzle exit. From the plot of 
variations in L presented in Figure 9 it is evident that the con- 
stant latent heat assumption in the basic equations is a reasonable 
assumption and any variation in L is negligible when compared 


to the effect of an equivalent variation in c. 


4, Nozzle Geometry 

The condensation process is a function of time, so it is ex- 
pected that the rate of expansion should have some effect on the 
"condensation shock. " Theoretical calculations were performed 
for the Nitrogen (I) expansion in a two dimensional nozzle for which 
the exit half-angle was varied from 9. 75 to 45 degrees. Nozzle 
lengths were selected so that a Mach 8 nozzle was used for each 
run. The results of varying the rate of expansion for Nitrogen (I) 


are presented in Figure 10 and Figure 11. In Figure 10 it is noted 





q 


that the faster expansion causes the vapor to supersaturate toa 
lower temperature before the onset of condensation. The "сопдеп- 
sation shock" is thinner for the fast expansion, it is (in physical 


dimensions) 3 cm for the 9 = 45° expansion versus 5 cm for 


1/2 


the 0 1/2” 9. 75° expansion. 6 is the nozzle exit half-angle. 


1/2 
This is somewhat misleading because, for Figure 11, in which 
all nozzles are referenced to a non-dimensional length for a Mach 
8 nozzle, the faster expansion tends to smear out over a broader 
relative length. The plot of the percent of condensed vapor in 
Figure 11 indicates that the higher the rate of expansion the less 
moisture is formed in the nozzle. It appears that for sufficiently 
high rates of expansion a condensation free flow can be obtained. 
This has been borne out by "condensation free" experiments in 


the past. 


D. Application to Metal Vapors 


Аз stated in Section V, the theory as developed herein should apply 
to metal vapors. The primary interest here is in the finding of a metal 
vapor that will condense in a highly expanded flow and that is suitable for 
laboratory experimental use. Hill (11) has predicted that sodium and 
potassium will condense; however, these vapors introduce critical han- 
dling problems. Copper, zinc, tin and lead appear to have possible experi- 


mental use, with copper being the preferred vapor of the four. 


Theoretical calculations were performed for a two-dimensional nozzle 
with a constant slope diffuser for the metal vapor input data listed in 


Appendix C. 





Ag. 


1. Copper 

Fortunately, copper vapor is theoretically predicted to condense 
for a fairly rapid expansion in a two-dimensional nozzle with a 20 
degree diffuser half-angle. The results are presented in Figure 12 
for an expansion in a Mach 10 nozzle from an initial pressure of 
approximately one atmosphere and an initial temperature of 3000°K. 
The "condensation shock" occurs very near the throat, and the change 
is very rapid with the "shock" spanning less than a centimeter of 
nozzle length. The average critical drop radius is approximately 
DES йй: ° cm and the maximum droplet formation rate is of the order 
of D drops per centimeter cubed per second. Again the onset of 
the condensation appears to be marked by a value of J = DA or 
w drops per ceritimeter cubed per second. Copper could be a 
Suitable choice for experimental work in rapid expansions because 
the large fraction of condensate formed in this theoretical solution 
indicates that higher rates of expansion and higher initial tempera- 
tures and pressures probably could be used without resulting ina 


"condensation free" expansion. 


2. Zinc 
Two computer runs were made for zinc, and neither run pre- 

dicted any detectable amount of condensate. Zinc (I) was computed 
for a 45 degree diffuser half-angle and initial pressure and tempera- 
ture of approximately one atmosphere and 1500°K respectively. The 
results of this calculation are presented in Figure 13 where a suffi- 
cient number of the calculated points are plotted along the isentrope 
to show that the expansion is essentially isentropic. This indicates 
a negligible amount of condensation. The tabulated calculations from 


the computer give a maximum g of 1. 48 x 10 at the nozzle exit. 
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The maximum nucleation rate obtained is J = 2 x 1012 drops per 
centimeter cubed per second and the critical drop radius is 

2х 10° cm. It was anticipated that a slower expansion witha 
lower starting temperature would predict a detectable amount of 
condensate. Zinc (IM) was run on the computer for a 20 degree 
diffuser and a pressure and temperature of one atmosphere and 
1300°K respectively. The results in Figure 13 again plot along 

the isentrope and indicate no condensate. The maximum g of 

4.8 x 1072 for the Zinc (II) run indicates some movement toward 
condensation. The maximum nucleation rate of J = 3 x dg. drops 
per cubic centimeter per second is obtained at a vapor temperature 
of 315? K and the critical drop radius is again 2 x 108. cm. The 

J = D value with no condensate formed lends support to the 

J = o criterion for the onset of condensation. A slower expan- 
sion of zinc probably would produce a condensation shock; however, 
this would be cut of the desired rapid expansion regime. From 


the results of the Zinc (I) and Zinc (II) runs it appears that for 


rapidly expanding flows zinc vapor may not condense. 


3. Tin and Lead 

Unfortunately, the computer runs for tin and lead failed to 
converge in the CFLOWS Subrcutine. The reason for this is not 
immediately apparent, and a further investigation will have to be 
made before any definite conclusions can be drawn. However, the 
few points that were calculated show promise that both tin and lead 


will condense in a rapid expansion. 





VIII. CONCLUSIONS 


An adaptation of the digital computer to the analysis of the condensation 


of a flowing vapor is possible, and the program for a pure vapor is pre- 


sented herein. Theoretical predictions of the condensation of nitrogen, 


copper and zinc indicate the following: 


I. 


The nucleation theory of Frenkel (13) provides a reasonable 
approximation to the rate of formation of condensate in a flow- 
ing vapor. 

The theoretical prediction for the condensation of nitrogen 
compares surprisingly well with the experimental results of 
Willmarth (3). 

A small drop size is the general rule. The critical drop radius 
is of the order of 10° cm and droplets of maximum growth 
rarely exceed ds cm. 

The nucleation rate controls the onset of condensation, and a 
prediction procedure based on this parameter might prove more 
usefu! than the one proposed herein. 

The growth equation plays a predominant role once condensation 
has started and a more accurate expression would be beneficial. 
The surface tension is the most influential parameter. Generally 
the value of o is uncertain and this leads to the largest source 
of error in the program. 

Values of specific heat and latent heat of vaporization are not 
critical, and treating these parameters as constants is a rea- 
sonable assumption. 

The condensation process is very sensitive to rates of expan- 
sion, and for rapidly expanding vapors a "condensation free" 


expansion is possible. 
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Copper vapor will condense in rapidly expanding flows and 
appears Suitable for experimental work in the condensation 

of metal vapors. 

Zinc vapor does not condense for high rates of expansion from 
one atmosphere pressure, and runs for lower pressures and 


temperatures are recommended. 
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Figure 1. Theoretical Results for Nitrogen Compared 
to Experimental Results of Willmarth (3). 
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Figure 2. Theoretical Prediction of the Condensation 
of Nitrogen (II) Compared to Experimental Results. 
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Figure 3. Typical Expansion of Nitrogen (I) Presenting Variation of 
Nucleation Rate, Maximum Droplet Radius and Critical Droplet Radius 
Along Nozzle Axis. 


Drop Radius, cm 





g, % 


-58- 


ТЕШЕ 
ТЕ Ж ЕТТТ 
Г T SN 

TTAB AO 

ARE A 
ИС 
АЕ 
LLL T T 
“ыр 
ТН LL L-I" 


ц 
КЕГЕ ЕТ 

ГУТ ЕЕ 
L TI I TL 
Г ГИ ТЕЕ 
mismas sama 
| |Z TO RP 
| TL TTE 
LIL O ITL L 
LLL I ЕТІС 


10 24 


x, Cm 

























Figure 4. Theoretical Prediction of the Effect of a Variation in the 
Accommodation Coefficient on the Condensation of Nitrogen (I). 





р/р. x 10^ 


-59- 










entirope 


8 10 12 14 16 18 20 
x, cm 


Figure 5. Effect of Choice of e as a Predictor of the Onset of 
Condensation for Nitrogen (T). 
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Figure 6. Effect of + 10% Variation in Surface Tension on the Theroetical 
Prediction of the Condensation of Nitrogen. 





-61- 


E... (6 A 

ше L 

— | £ el Б 

— 8H 
‚10, 


"HEIN 
! ! 





2 


p, dynes/cm 





Figure 7. Effect of Variations in Surface Tension on the Predicted 
Super saturation of Nitrogen. 
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Figure 8. Comparison of 10%. 1d 50% Changes in Specific Heat to a 10% 


Change in Surface 'l'ension for Nitrogen (I). 
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Figure 9. Comparison of 10% ana 50% Changes in Latent Heat of Vaporization 
to a 10% Change in Surface Tension for Nitrogen. 
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Figure 10. Theoretical Expansion of Nitrogen (I) for a Two-Dimensional 
Nozzle with Varying Exit Half-Angle, 9 1/2 
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Figure 11. Theoretical Expansion of Nitrogen (1) in Mach = 
8 Nozzle of Lengths 35, 16.5, and 6 cm. 
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Figure 12. Theoretical Prediction of Condensation for the Expansion 
of Copper Vapor ina Two-Dimensional Nozzle. 
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Figure 13. Theoretical Prediction of the Expansion of Zinc (I) 
and Zinc (II) Vapor Indicating No Condensation. 





Аррепдіх А 
DERIVATION OF THE NUCLEATION RATE EQUATION 


This is in no way intended to be a complete derivation of the nuclea- 
tion rate equation. The full derivation can be found in complete detail 
on pages 366 through 400 of Frenkel (13). A scant derivation is presented 
in Stever (12), pages 532 to 540. 


There appears to be an error of 7 in the final equation presented 
by Frenkel, so enough of the derivation will be reproduced here to clear 
up this discrepancy. The solution is a statistical problem of kinetic 


theory. 


In an equilibrium distribution of droplets where $ < ШЕ the fol- 


lowing relation must hold. 


N S о = М 8 (А.1) 


noq ED n - 1 


where ра апа ali are the thermodynamic potentials of the vapor and 


liquid drops, respectively. N. and N are the number of droplets 


n-1 
containing n and (n - 1) molecules, respectively. S and S ME 
the surface areas of droplets containing n and (n - 1) molecules, respec- 
tively. o is the probability that one molecule leaves a unit surface 


area of a droplet containing n molecules. 


В, Ж is the probability that one molecule condenses on a unit surface 


area. Assume о, апа B to remain valid for the non-equilibrium case 


Piia 
tion f The equation will no longer balance so (A. 1) become s 


1 
«d Фф.) and modify the distribution N. to a non-equilibrium distribu- 
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Е О З D NEN IU (A. 2) 


where J represents the increase per unit time of the class (n) due to 
condensation of the vapor on the surface of droplets in class (n - 1) over 
the number, which, owing to evaporation, pass from the class (n) to 


the class (n - 1). 


The rate of change of the number of drops of a given class is deter- 


mined from (A. 1) and (A. 2) to be 


= INS (А. 3) 


Frenkel employs the method of Zeldovich which utilizes the feature of 
the sharp maximum of Ad at n* to solve the equation. By considering 
only values of n > 10 say, quantities appearing in the equation will 

only vary slightly when n is changed by 1 and can be treated as functions 
of a continuous variable. Equation (A. 3) can then be integrated and the 


result is presented as equation (27b) on page 396 of Frenkel (13). 





tal AS (n*) Y 
= * = 
J М, Dín*)| exp ET | 2 ET (А. 4) 
where N = Бе. number of molecules in the supersaturated vapor 
a Sa 1/3 2/3. 2/3 2/3 
D(n*) = S(n*) 6 = diffusion coefficient = (4r ) 3 MS n* B 

p EN 32 
Vv 27 mkT 
Аф = =" С re? 


el ж 
y = gs ф Ип 


substituting these values into (A. 4) produces the equation 


2 7 2 1 (Ф нд = Фф). 


ror* 
J = 2г*° |р. je елес (А. 5) 


SKE | V3 mn* 
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or, writing in the form of Frenkel's final equation (28a) 


n*J ope? Ше _ 47 E 1 Gig j | 


М КТ 
n 





= ann 
3kT 3 m 2 (A. 6) 


In comparing (A. 6) with Frenkel's equation (28a), it is noted that Frenkel 
does have an extra factor of 7 in the ccefficient as noted by Courtney (9) 


and others. 


The desired equation for this analysis is obtained by substituting 























20 V.. 
the value of (Siig - Фф.) from the critical radius equation r* еч 
11д у 
and the value n* = = 27 ros into equation (A. 5) 
liq 
2 
2| p 
= * === eo 
J = 2r m| [EXP 
s lle: c 
орж liq | KT ЗЕТ тіп 
2 2 
NE 4/39. exp - S70 r* 
ar кт Ула Ута er (А. 7) 
To eliminate the necessity of calculating molecular values of 
Vu m and m = E , the equation for computing J in this 
E AT, A 
analysis is written 
KT) pp (TN, P ` 3kT | 











Appendix B 


IBM 7090 DIGITAL COMPUTER PROGRAM FOR 
CONDENSATION IN HIGHLY EXPANDED FLOWS 


Written in Michigan Algorithmic Decoder (MAD) Language 


This program is valid for a pure vapor with constant specific heat 


and constant latent heat of vaporization. 


The cgs system of units is utilized throughout and cgs values for 
the universal gas constant, Boltzmann's gas constant, and Avogadro's 


number have been included within the program. 


The symbolic names of parameters used in the program are listed 


and defined at the beginning of each applicable section of the program. 


The units for and input data required for utilization of this program 


are presented in Appendix С. 


ыт 





Ес 


SCOMPILE MADsPRINT OBJECT »PUNCH OBJECT MAIN QOLMAIN 000 


START 


AERONAJTICAL ENGINEERING 


е (44-5 2-6 E < - w “m. Oe ` — =... or æ vame o 








R NOZZLE CONDENSATION RESEARCH 
R 
— 7 = MAIN PROGRAM. I тала аса 
R 
RRESEARCH-- Je GRIFFIN PROGRAMMING-- Le HARDING JULY 1963 
R 
PROGRAM COMMON GAMMA» MUs Bs Cs CPs Ls RHOL> SIGMA . 
R 
Eu ATA. =o  — "EE S за rs Q as 
R ЖҰ АНТЫН ee 
RTHIS DATA SHOULD INCLUDE Y Ba 
R 
R PZERO = INITIAL PRESSURE | малта ы 
R RHZERO = INITIAL DENSITY "NN dE 
R TZERO - INITIAL TEMPERATURE - | m 
R PE - ТБЕМПКОРТС ЕХ PRESSURE su US ЖЕ БЕРСӘ NN 
R INITIAL APPROXIMATION TO THE SATURATION PRESSURE. 
R 
TR GAMMA = RATIO OF THE SPECIFIC HEATS "0110 2 
R MU = MOLECULAR WEIGHT OF THE VAPOR m ии 
СЕ СР = SPECIFIC HEAT AT CONSTANT PRESSURE 
R L = LATENT HEAT (NOTE VAPOR PARAMETERS SUBROUTINE) I 
R SIGMA = SURFACE TENSION (NOTE VAPOR PARAMETERS SUBROUTINE) 
R 
ЕС 7 - SATURATION CURVE CIN TANI no — LLLI 
R B = SATURATION CURVE CONSTANT E 7и 
R ALPHA = ACCOMODATION COEFFICIENT 
R 
^R DELX = INCREMENTAL STEP IN Xe THE NOZZLE CONDITIONS WILL. 
R _ BE COMPUTED AT XCON+DELX»XCON+2*DELX>... WHERE 
R XCON IS THE CONDENSATION POINTS 
R XRANGE = LENGTH OF INTERVAL» STARTING AT THE CONDENSATION 
R > POINT» IN WHICH THE FLOW CONDITIONS WILL BE 
R COMPUTED. | _ Е 
^R XPOINT = THE CONDENSATION POINT IS ASSIGNED A SUBSCRIPT OF 
R ZERO» ONLY THOSE POINTS WITH SUBSCRIPTS XPOINT» | 
R 7 7 77 Z2*XPOINTseee WILL BE PRINTED. 
R 
R TRANGE = PROGRAM WILL LOOK FOR A CONDENSATION POINT IN THE 
R INTERVAL (TSATsTSAT~TRANGE)» WHERE TSAT IS THE | 


С ^^ COMPUTED SATURATION TEMPERATURE. 
R DELTAT = INITIAL STEP SIZE USED IN LOOKING FOR CONDENSATION 
СЕ 


“POINT IN THE INTERVAL (TSATsTSAT-TRANGE) 
R PERCNT = IF EPSLON(T) EXCEEDS PERCNT IT IS ASSUMED THAT 


u. 


—— — — M—————— + 








R "CONDENSATION HÁS STARTED. 
R 
~ PRINT COMMENT $1 NOZ 
1ZLE CONDENSATION RESEARCH %- x ЕЕ 
PRINT COMMENT $O | INITIAL CONDITIONS$ 
PRINT RESULTS PZERO»RHZERO»TZERO»PE 
— "PRINT COMMENT $0 ~ “CONSTANTS DEPENDING ON THE VAPORS 
PRINT RESULTS MUsCP»GAMMA»SIGMA »L 
"PRINT COMMENT 30 °° SATURATION CURVE CONSTANTS$ 
PRINT RESULTS C,B 
— —PRINT COMMENT 50 — — PROGRAM PARAMETERSS$ A ^ "^" 7^ 


PRINT RESULTS DELX»XRANGE »XPOINT»ALPHA 
~ EXECUTE NOZZLE e(STHROATSs00esASTAR) Dp 5 





52 


$3 


sod 


FIND SATURATION POINT BY MATCHING ISENTROPIC PRESSURE 
FORST TON TU TRES AMOUR O NS UR VEA ES SURE 


A 4x 


pspr 

2-САММА/(САММА-1») 

== SECOS. 1P)) 

Bear Ze ROe( CTZTZERO) e ) 

TEST=+ABS. 11. -PSAT/P)=-+.001 

MeN VES TCSTeleGes | RANSEER 10752 

F=P+,5#(PSAT-P) 

IESSNSEEX LO 54 

TSAT=T ; 
EXECUTE IFLOWSs(ASTAR»CP»GAMMA,PZERO»RHZERO»TZERO) 
EXECUTE IFLOWes (ASAT»MSAT»PSAT»RHOSAT»TSAT»USAT) 
Z=MSAT 

EXECUTE NOZZLE«(SINVERSSsZsASAT ) 

WHENEVER Z(1)»*E«* 14 » TRANSFER TO 8S3 


PRINT COMMENT 34 THE SATURATION POINT FOUND IS NOT INSIDE 
] THE NOZZLE» THE PROGRAM WILL CONTINUE+5 
XSAT=Z 
MDOT=ASAT*RHOSAT*USAT 
PRINT COMMENT >2 COMPUTED SATURATION CONDITIONSS$ 


PRINT RESULTS PSAT»RHOSAT »sTSAT 
PRINT RESULTS MSAT»USAT ~ 
PRINT RESULTS XSAT»ASAT 

PRINT RESULTS МОСТ 


FIND CONDENSATION POINT BY T BY MATCHING ISENTROPIC FLOW 
SOLUTION TC THE NUCLEATION THEORY SOLUTION GIVING THe 
AMOUNT OF CONDENSATE AS A FUNCTION OF AN ASSUMED FLOW- 


AD AN AN 


DELT=DELTAT 
ШЕРЕЗ0» 
TMIN=TSAT-TRANGE 
WHENFVFR Tell Nie lei Oe » TEIN=08 
— "PRINT COMMENT $2 PARAMETERS DETERMINING СОМОЕМ5АТІО 
1N POINT SEARCHZ VPN mu 
PRINT RESULTS TSATsTMINsDELT 
PRINT RESULTS PERCNT 
“ 
R 
ПЕТА ТУ I 
[ЕЕ =О he. Е КҮРЕ. 
^ T*5T-DELT 
WHENEVER T eLe TMIN 
— 2 ass —— 
ITERF=ITERF+1. | 
“рТ к” a 
WHENEVER ITERFeLelOe»TRAN>FER гу 54 BER 
WHENEVER Ze e00esTRANSFER ТО 57 
PRINT COMMENT $0 TEN PASSES OVER TEMPERATURE RANGE SPECIFIED» 
i NO CONDENSATION POINT INDICATEO.$ 
TRANSFER TO START | 
o END OF CONDITIONAL 
EXECUTE IFLOWe (AsMsPsRHOsT 5U) 
Z=M 
EXECUTE NOZZLE+(SINVERDDILOIA) 
WHENEVER Z(1) eEe Ov 
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EX. 
| TRANSFER TO $6 nn ШЕ m 
Еу PRINT COMMENT $0 TEN PASSES OVER THE FULL RANGE OF THE NOZZLE 


_ls NO CONDENSATION POINT INDICATEDeS | 
"TRANSFER TO START IT C C d 
END OF CONDITIONAL 

кос с аена 
EXECUTE ОАВ з Ееее рава 2) _ 

Т БЕСА =А 

В 

ИНЕ Г „(Ту 
TDROP=C/ (B-ELOG+(P)) 
RHOL=FRHOL+ (TDROP) ` 

 SIGMA=FSIGMAsLTDROP) — _ _ = С . 
EXECUTE NUCLEle (AsDELXsMDOT sP oT sRADsNDOT » DELG) 

R 

R COMPUTATION OF EPSILON (T) 

R 

- EPSLON=(GAMMA=(1+/(M*M)))/(GAMMA-10+) 
__WHENEVER MeLele»EPSLON=le 
EPSLONSA*(1.-(L/(CP*T))*EPSLON) — 
EP SLON=EPSLON*(DELG/DELA) 
_ EPSLÖN=.ABS. EPSLON ` 

R 

~ WHENEVER EPSLONeLePERCNTs TRANSFER TO 55 p 
T=T+DELT 

O DELTA, 5DE nn m созт 
ITER=ITER+1, 

WHENEVER ITEReLeLITER»sTRANSFER TO S5 


"w s“ >. . w. аа w r" NEUES e$ 5 9 + + + == Чие e re > Ve ры 


-— o o 


a 


un ree NT TUR ++ + . a + 





SS E E м E 2) 
|Á— —— — —. — — — -— .- —— — —S>——ə—— 


 — [nF E — - ———,.— oo -— -- — -e — -- — _ ~ - o 


LITER HALF-INTERVAL STEPS HAVE. BEEN PERFOR 
A CONDENSATION TEMPERATURE HAS BEEN FOUND.» 


— 


AND 


m mm 


PRINT COMMENT $0 HAVE FOUND TEMPERATURE THAT SATISFIES THE CO 


en а 


INDITIONS FOR THE CONDENSATION POINTS 
PRINT RESULTS T»DELT 


TCON=T 
EXECUTE IFLOW. (ACON»MCON »PCONSRHOCON»TCON»UCON) ` 


¿ MGON — 
EXECUTE NOZZLE.(SINVERSS»Z,ACON) 


ИШЕМЕУЁК ZOOL) EsTe + TRANSFER TO S8 


= - -- - — —v ee nn un _— sm- et ee o e — 


PRINT COMMENT %4 THE CONDENSATION POINT FOUND IS NOT INSI 

1DE THE NOZZLE» THE PROGRAM WILL CONTINUEeS | 
58 XCON=Z 

Er Je M 


TDROP=C/(B-ELOG+ (PCON)) 
—RHOL=FRHOL+(TDROP) 7 
SIGMA=FSIGMA+(TDROP) 

EXECUTE NUCLE1le (ACONsDELX,»MDOT»PCÖN»TCON»RAD»NDOT sDELG) 
G=DELG 
"PRINT COMMENT $4 777 | ‘COMPUTED CONDENSATION CONDITIONSS 
PRINT ee Ыса. 


— cm. m amom a oT amm a 


PRINT ur: XCON» eon 
"PRINT RESULTS SIGMA,RHÖL»L»CP ` 

PRINT COMMENT $ NUCLEATION THEORY QUANTITIESS 
PRINT RESULTS RAD,NDOT sDELG»G 


R 
RINITIALIZE FOR NOZZLE COMPUTATIONS STARTING AT THE CONDENS= _ 


- = Se 





R DELA 


- PPOINT=0 


^7 DELU-.0001*0CON 
СЕЕ ОЕ 


^». U(1)sU-*DELU | 


ne 


RATION PCINTs THE VECTOR X CONTAINS THE POINTS ALONG THE 
Rolo OF THE NOZZ2LE AT ЕТЕН THE QUANTITIES DESTRED SARE wie 
RbE COMPUTEDe NOTE THAT IT IS NOT NECESSARY THAT THE SPACING 
RBE UNIFORM SINCE A VECTOR DELX IS PROVIDEDe AT EACH POINT 
ROF THE MESH WE WANT TO CALCULATE 


R P = VAPOR PRESSURE 

R RHO = MIXTURE DENSITY 

R T = TEMPERATURE 

Q J = FLOW VELOCITY 

R M = MACH NUMBER 

R DELG = CONDENSATE INCREA -E SINCE THE Al) Meche cen 
R о - ЕКАСТІОМ» БҮ МА55» ОҒ THE VAPOR THAT IS IN THE 
R CONDENSED >TATE | 
ROcLR(T) = THE INCREASE IN RADIUS» SINCE THE LAST POINT, CF 
R THE DROPS INITIALLY FORMED AT POINT I 

R RADCIY = "RAD USsOr DROPS SUNT TT ALLY seOR MED sae 100) 
RNDOT(I) = THE NUMBER OF DROPS PER SECOND THAT ARE NEWLY 
R FORMED FROM X(I-1) TO XI) 


R SURFACE TENSION 

R EOOD DENSITY 

R SPECIFIC HEAT AT CONSTANT PRESSURE 
R LATENT HEAT 

RALSO THE IDENTIFYING GUANTITIES 

R 

б 

R 


ж 
E 
O 
i 
" u.m Hu 


X POINT OF EVALUATION 
DELX = DISTANCE FROM LAST MESH POINT 

A NOZZLE AREA | | 
AREA INCREASE FROM LAST MESH POINT 
RONLY THE POINTS» THEIR INCREMENTS AND THE INTEGER INDEX М 


-— 


RARE KEPT IN CORE, ALL THE REST OF THE INFORMATION IS MAIN- 


RTAINED ON TAPE« ALL QUANTITIES FOR POINT Ns HOWEVER» ARE ` 
RREQUIRED TO OBTAIN THOSE FOR POINT N*l»s HENCE THEY ARE 
 RTEMPORARILY IN CORE AS UNINDEXED QUANTITIES WITH THE ABOVE. 
 RNAMES. DM е Ф 

R 


INITIALIZE TO CONDENSATION POINT OE S D ы 
R 


ыа э 
P=PCON 
 RHOZRHOCON 
U=UCON o 
 M=MCON | ili 
X=XCON 
А 

EXECUTE NOZZLE+ (SAREAS >) XCON-DELX>A(1)) 





= a = 


^ DELAsA-A(1) 


N=0 


EXeGUTE TOCTRL« 


R 


—— ——— €———————— 


WHENEVER U( 1) eLeOe 


~~ PRINT COMMENT $4 ~ THE INITIAL APPROXIMATION TO THE VELO 


ICITY AT THE NEXT MESH POINT IS МЕСАТІУЕ»% 


TRANSFER TO S34 
END OF CONDITIONAL 


—R————-- 


— ct ——. --- --- - .------ nn —o Ée; ——— —— — — —— ——  — ——— — 


—— A A — —————— — 2. ....- - — —ÓM —— -- 


A Ha A - --- ------“..-- .— 


———————————— аа . =. —— 
с — PX T K | 


% 





U^ 
м 
к 
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R GEGIN ITERATION -ON DELG UNTIL NUCLEATION THEORY AND 
R FLOW EQUATIONS AGREE. A MAXIMUM UF GITER ITERATIUNS ARE 
R ALLOWED PER MESH POINT. 


Meet =o ТЕЕ 


N=N+] 
WHENEVER NeGelCQ0 .... 
PRINT COMMENT $4 NUMBER OF MESH POINTS LXCFEDS THE PRO 2 
1GRAVMED LIMIT.4$ е ° 
TRANSFER TO S34 ° 
END OF CONDITIONAL m .... 


Pee xX Ch) = 06 ix 
X(N) =X(N-1])+DELX(N) 


w S EVER XG XLIMIT Фоо е 
PRINT CCMMENT $4 THE NEXT POINT WOULD Reaper ЕЕ Š 
ІІ БЕБТЕЕ 5% : 
TRANSFER TO S34 . 
END EOE ONDI IONAL 2-7 
EXECUTE МО271Е.(ФАЗЕА>ЭХ(М)»А(1)) 
WHENEVER A Om Фоо о 
PRINT COMMENT 34 NOZZLE SUBROUTINE ERRUR INDICATIGN - Я 
BEESHOESPCINTODOES NOT LIE INSIDE NOZZER.S e 
Po ieee SULT Sex CN) ° 
DRENSPER 70 5238-2-77 УМ” кке ° 
END GF CONDITIONAL . 25822 


MEA SA 
CELG(1)=DELÓ ' 
G(1)=G+DELG(1) Е 
EXECUTE CFLOWe (A(1)»DELA(1)>G(1)»DELG(1)>»P>»RHU»>T>U9M) 
L=FLe(T(1)) 
TDROP=C/(B-FLOGe(P(1))) 
RHOL=FRHOL + (TOROP) 
SIGMA=FSIGMA.(TDROP) 
TEST=G(1) 
CEXECUTE NUCLE2Z> (X(N) sDELXIN) sACl] J sDELA( 1) 3P(1 Is Reet ae 
1U(1) 9M 1) sMDOT sALPHA sN sRAD »DRAD sNDOT sDELG(1)) 
~ G(1)=G+DELG(1) 
TEST=.A55.(1E51-6(19) 
WHENEVER TESTeLe «00001 s TRANSFER ТО 532 
Pier alr R= ls 
~ WHENEVER ITEReGEeOesTRANSFER TO 531 
PRINT COMMENT % 2.7272 9792 9 9 9 D 
Ex coc IU 
DEIMXEDELX(N) 
Esca D 
DELASOELADI) 
—P=P(1) 
RHO=RHC (1) 
UE S 
DELU=SU(TI-U 
9-0(1) 
M=M(1) 
BOEMG=DFLG( 1) 
G=G(1) 
THROUGH S339FOR I=0>l1l>Il.Go[N-1) 
RADCI)SRADCI)TDRAC(I) 


e A -A»M- +` аға. ш» ч? у + EIU -. 


RÁ 


NECE остер. 
R 





DD DAD 
M 


X PRINT RESULTS RAD(O)e« RAD(N-1) 


" PPOINTsN*XPOINT 


INTEGER 1I+N+PROINT | 


| an 
TRANSFER TO 529 
R 


TRANSFER TO START 


INTERNAL FUNCTION FOR OUTPUT OF THE FLOW QUANTITIES AT 
MFSH POINT Xe LATER THIS SECTION MAY ALSO SAVE THIS 
INFORMATION FOR FURTHER PROCESSING. 


INTERNAL FUNCTION 
R 
_ ENTRY TO IOCTRL» 
R 
WHENEVER  N«E,PPOINT. e | | | 
PRINT 2-1 54 NOZZLE CONDITIONS -- $ 
PRINT RESULTS N»XsDELX»AsDELASP SRHOS T»USMSRADIN) SNDOT (N) » 
1DELG 5G 
WHENEVER NeGe0 


END OF CONDITIONAL _ 


w w —ww s a a ` - - 


END OF CONDITIONAL 


= аа ааа - —- — 2... 322022. - = 


FUNCTION RETURN 


END OF FUNCTION 
R 


DIMENSION DELX(1000) »DRAD(1000) sNDOT(1000) »sRAD( 1000) »X(2000) 
DIMENSION A(1)»DELA(1)»DELG(1)9G(1)»M(1)»P(1)»RHO(1)>T(1)» 
Meee yee) 
VECTOR VALUES XPOINT = 
VECTOR VALUES LITER = 
VECTOR VALUES GITER = 
R 


— -— — - . — .. - - --. — 


END OF PROGRAM 
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ILE MADsPRINT OBJECTsPUNCH OBJECT IFLOWOOLIFLOWOO* 
RS AERONAUTICAL ENGINEERING EN. sss 
R “NOZZLE CONDENSATION RESEARCH 
R 
R I ISENTROPIC FLOW Нак 
R 
Te EARCH-- Je GRIFFIN PROGRAMMING-- Le HARDING JULY 1963 


‚КО торт FLOW EQUATIONS FOR COMPUTING PRESSURE, MACH NUMBER 

RAREA» VELOCITY AND DENSITY GIVEN THE TEMPERATURE AND THE 

RFOLLOWING PARAMETERS 

ASTAR - THROAT AREA OF NOZZLE 

“ССР = SPECIFIC HEAT AT CONSTANT PRESSURE 

GAMMA = VAPOR EXPONENT 

2277 Р2ЕБКО - INITIAL PRESSURE 

RHZERO - INITIAL DENSITY 

"— TZ2ERO - INITIAL TEMPERATURE = 


- + @% . a. 6 -— - en —.. vw A a A AX EEE n — s 


DAA AA AA 


TeEX PERNA EUNGT ION ТАТ,А2УАЗУДАДУАТАС-) 9 = 
07 “THIS ENTRY PROVIDES FOR PICKING UP THE PARAMETERS LISTED 
ABOVE» THEY ARE THEN SAVED INTERNALLY FOR LATER USE. 


ENTRY TO IFLOWS. 

оолу у eae 
СР-А2 
GAMMA=A3 
PZERO=A4 
лра н ccc cc j ——єу. 
TZERO=A6 
FUNCTION RETURN ЕУ I 


A 


-yrr -n me чыз, “ча - т, ` ~ e > . m. # ~ 2.28.4 - aw < — w umasa ADE a pp —— ee + rm s е 


THIS ENTRY COMPUTES THE UNKNOWN QUANTITIES 
А1 NOZZLE AREA 
A2 MACH NUMBER 
АЗ PRESSURE 
Aá = DENSEPY I ( PiO en ne 
A6 VELOCITY 

IN TERMS OF THE PARAMETERS SAVED FROM THE LAST CALL 

OF 'IFLOWS*? AND 
PAS = TEMPERATURE 5 БЕСЫ m 

NOTE THAT THE LIST OF ARGUMENTS IS ALPHABETICAL. 


л) 


| 


jn m n wn n 


д) д д] л) л) л л Ds 


ENTRY TO IFLOWe 
UD o 
R PRESSURE 
~ EXP=GAMMA/ (GAMMA~1  ) 
L ZERO; Ро EXP 
 A3=PZERO*2 ALS 
R WEEDCITY AND MACH NUMBER. 
Z= 2 [Z RO-T) 2 
A6G=SQRT+. (CP*Z) 
-2:2/(САММА-1.) | 
AZE SORT: (Z/T) 
R DENSITY AND NOZZLE AREA " | Ten ^. uM 
ЕХР-1./(САММА-1») 
7-251%44%5ЖАОЖА?/ЕХРО-77 
AG=RHZERO/(Z «Pe. EXP) 
“-“ЕХРЕ.5%( САММАҒТІ» ) “ЕХР 


0 A m  - DO @ — 


wa е - - a. A oo тш AAA жы” лынған 
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.%2/(САММА-1.) 

ЖЕТЕ EXP) /AC 

A 
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Eu MPILE MAD,PRINT OBJECTIPUNCH OBJECT CFLOWOOICFLOWOO* 


К AERCNAUTICAL ENGINEEKINOG 

R NOZZLE CONDENSATION RESEARCH 

К 

К CONDENSATION FLOW 

R 

RRESEARCH-- Je GRIFFIN | PROGRAMMING-- Le HARDING JUL aso 
R 


RCOMPUTATION OF VALUES OF PRESSURE DENSITY s TEMPERATURE » MACH 
RNUMBER AND VFLCCITY THAT SATISFY THE FLOW EQUATIONS THAT 
RAPPLY AFTER THE ONSET OF CONDENSATIONe IT IS ASSUMED THAT WF 
RHAVE COMPUTLD THE CONDITICNS AT X AND DESIRE THOSE AT X+DELX 
RTHE ARGUMENTS REQUIRED ARE 

R DELA = AREA INCREASE FROM X TO X«DELX ` 


R A = AREA AT X+DELX 

R DFLG = CONDENSATE MASS INCREASE FROM X TO X+DELX 
R G = CONDENSATE MASS AT X4#DELX Т5 

R P = PRFSSURE AT X 

R | ВАС s DENSITY AT X Е 

R T = TEMPERATURE AT X 

R UE VELOCITY AT X 

R M = MACH NUMBER AT X 


RTHE PARAMETERS CP,L AND MU ARE REQUIRED FROM PROGRAM COMMON. 
pant FUNCTION FLe(T) MUST BE USED HOWEVER SINCE |) ERA ГО 
RT IS NECESSARYe THE UNKNOWN QUANTITIES ARE RETURNED IN P(1)» 
ReeOt lyst (lysUtl) AND М(З2 Il 1S ASSUMED THAT VU) ИИ ЕЕ 
CENTINITIAL APPROXIMATION TO THE VELOCITY Al (te NES 
REOLNT ON ENTRY 


R 
EXTERNAL FUNCTION (A»DELA»G»DELG»PA»RHOA»TAsUA»MA) 
R 
PROGRAM COMMON GAMMA» MU» B» С» СР» 1» RHOL» SIGMA 
R 
ENTRY TO CFLOW. | 
U=UA(1) 
DELU=U-UA 


MDGG=DELG/(1l.-6) 
DQA-DELA/A 

PERSIA) 
R 

DQUZDELU/U 
R 
R SOMPUTATION ОЕ ВЕСТЖАМО Т ЕСО ЫЕ ENERGY EQUATION 
R 
BEEISDIESDEPG-USDEDLU)ZCP 
T=TA+DFLT 

L=FLe(T) 


R COMPUTE DELRHO/RHO FROM CONTINUITY EQUATION 
— m 

DQRHO=DQA+DQU 

DQRHO=-DQRHO 


R 
—R — 7 СОМРОТЕ DELP/P FROM THE EQUATION OF STATE s 
R eee UE NE ms Pin ae 
—DaPzDORHO*(DELT/T)-MDGQGO ^ 
R 


R WE CAN NOW COMPUTE THE VELOCITY CORRESPONDING TO 
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R THESE DIFFERENCE QUOTIENTS FROM THE MOMENTUM EQUATION 


DELU=-(DQP*R*T*(1+-—G))/(MU*U) 

NEWU=UA+DELU 

ТЕ5Т-.А55.(1.-МЕМ(//1/) 

U=NEWU 

WHENEVER TESTeLeEPSLONs TRANSFER TO S3 
R 

R ITERATION LIMIT HAS BEEN REACHED WITHOUT CONVERGENCE 
R 
— PRINT COMMENT $  'CFLOW' — ITERATION FAILED TO CONVERGES 

EXECUTE ERROR, 

R 

R CONVERGENCE CRITERIA SATISFIED 

R 

TA(1)=T 

UA(1)=U 

DQRHO=1+-DQRHO 

DQP=1.-DQP 

PA(1)=PA/DQP 

RHOA (1) =RHOA/DQRHO 

NEWT=R#GAMMA#T/MU 

MA(1)=U/SQRTe (NEWT) 

FUNCTION RETURN 

R 

VECTOR VALUES R=8.314E+07 

VECTOR VALUES EPSLON=1.E-06 

INTEGER 1 

R 

END OF FUNCTION 
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BeeMPTLE MAD,PRINT CRJECT»PUNCH OBJECT NUGEDTOOCINUCETOG 
R AERONAUTICAL ENGINEERING 
R NOZZLE CONDENSATION RESEARCH 
R 
R NUCLEATION THEORY = 1 
R 
RRESEARCH-- Je GRIFFIN PRUGRAMMING== Le HARDING JULY Sc 
R 


RCOMPUTATION OF THE MASS DUE TO INITIAL CONDENSATION AS GIVEN 
RBY NUCLEATION THEORYe THE FOLLCWING QUANTITIES MUST bE AVAIL- 


RABLE. 
КІ!) CONSTANTS 
R NA = AVOGADROS NUMEER (MOLECULES/GRAM MCLE) 
R K = BOLTZMANN CONSTANT (DYNEXCM/DEGREEK) 
R KR = UNIVERSAL GAS CONSTANT (DYNE*CM/GRAM MOLE*CEGREEK) 
R(2) PARAMETERS 
R БОЕ САТОВА TON CURVER CONSTANT 
R C = SATURATION CURVE CONSTANT 
R DEX = DELTA X FOR VOLUME 
R MOOT = RACSAT#ASAT#USAT 
R MJ = MOLECULAR WEIGHT OF VAPOR 
R(3) VARIABLES 
К A = NOZZLE AREA 
R PZZSERESSURE 
R T = TEMPERATURE 
R RHCL = LIQUID DENSITY DEPENDENT ON TEMPERATURE 
R SIGMA = SURFACE TENSION DEPENDENT ON TEMPERATURE 
R 
EXTERNAL FUNCTION (A19A29A39A49A55A6 9A7 9AB8) 
R 
PROGRAM COMMON GAMMA» МО» &s Cs CP» Ls RHOL»s SIGMA 
R 


RCOMPUTE THE RADIUS OF THE DROPS CONDENSING UNDER THESE 
RCCNDITIONSs THEIR NUMBER AND HENCE THE CONDENSATE MASS 
RFORMED IN A VOLUME ELEMENT OF LENGTH DELXe THE PROGRAM 
RCOMMON PARAMETERS RHOL AND SIGMA ARE ASSUMED TO CORRESPOND 
RTO THE DROP TEMPERATURE FOR THE GIVEN FLOW CONDITIONS». 
R 

ENTRY TO NUCLEI. 

A A] 

DELX=A2 

ИВО r= A3 

Pp=A4 

[EE 
QU 

R COMPUTE RADIUS ( CMe ) OF THE DROPS THAT ARE CONDENSING 
R 

RADIUS2(2,*SIGMA*MU)/(RHOL*T) 

RADIUS=RADIUS/(R*(ELOGe(P)-B+C/T) ) 


COMPUTE THE NUMBER OF DROPS OF THIS SIZE THAT ARE 
CONDENSING PER CENTIMETER CUbED PER SECOND. 


TOAN 
1 


TEXP=RADIUS*RADIUS/K 
77 ПЕХБЕ-(4.1887902%56ІОМАЯТЕАР ЛТ ТТ Тат ен 
N=SORTe(SIGMA#MU/NA(1)) 
ШШ РОО ЛҮҮ у 77777 
N=P*N*EXP.(TEXP)/(K*T*RHOL) 
"^ Nz.7978846*N | 


sa- o mvae ouu aamu cmu e ew tee wk umm m eo o mn 4 (009 9 A o5 ey 
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NDOT =N*A*DELX 
R 
R COMPUTATION OF PERCENT OF LIQUID MASS» DELTAG 
R 
DELG=4+.1887902*RHOL*NDOT*RADIUS 
DELG=DELG*RADIUS*RADIUS/MDOT 
EREMO 
_AS=RADIUS 
AT=NDOT 
AB=DELG 


R 
VECTOR VALUES МА-6.027Е%23»6.027Е%03 
VECTOR VALUES K=1.379E-16»1.379E-06 


d 
Ooo O e e шш er rn a 


wars 
` ——— 2... ee аршалы DU ee е 


МЕСПОКТУЛМЕЦЕС  К58%321 6501 
К 
END OF FUNCTION 
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SCOMPILE MAD+PRINT OBJECT+PUNCH OBJECT | NUCL2001NUCL200* 





R — LLAERCNAUTICAL ENGINEERING — e 
R NOZZLE CONDENSATION RESEARCH 

R 

R |.  — NUCLEATION THEORY - 2 

R 

RRESEARCH-- Je GRIFFIN PROGRAMMING-- Le HARDING ДЕ ЖЕТСЕ 

R 


—RCOMPUTATION OF THE QUANTITIES ARISING DUE TO CONDENSATION 
RUNDER THE FOLLOWING CONDITIONS 





R X = NOZZLE COORDINATE 

К DELX = INCREMENT SINCE LAST MESH POINT 

R A = NOZZLE AREA AT X 

R — DELA - NOZZLE AREA INCREASE SINCE X-DELX = 
R P = PRESSURE AT X 

R RHO = MIXTURE DENSITY AT X 

R T = TEMPERATURE AT X 

R Ul = VELOCITY OF FLOW AT X 2 
R М = MACH NUMBER AT X 

R___MDOT = TOTAL MASS FLOW RATE i 
R ALPHA = ACCOMODATION COEFFICIENT 

R N = MESH POINT NUMBER SINCE CONDENSATION» THE POINT 

R OF CONDENSATION IS MESH POINT ZERO. 

RTHE QUANTITIES TO BE COMPUTED ARE 


i RADIN) RADIUS OF NEW DROPS FORMED AT Xe RAD(I) IS RADIUS | 
САТ. ASDELX Or DROPS SIND ОБ AT MESH POINT 
I AND IS NECESSARY INPUT ) 


жала.» x. 





RNDOT(N) = NUMBER OF DROPS OF RADIUS RAD(N) FORMED BETWEEN _ 
R 2 X=DELX AND Xe ( NDOT(I) IS THE NUMBER OF DROPS 

R INITIALLY FORMED BETWEEN MESH POINTS I-1 AND I AND 
КЕ JS NECESSARY INPI 

RDELR(O) = DELR(I) IS THE RADIUS INCREASE BETWEEN X-DELX AND X 
R THRU OF THOSE DROPS INITIALLY FORMED BETWEEN MESH 
RDELR(N-1) POINTS I-1 AND I» THESE QUANTITIES MUST ALL BE 

R COMPUTED. ) 

R DELG = CONDENSATE MASS INCREASE FROM X-DELX TO X DUE TO 
R NEWLY FORMED DROPS AND GROWTH OF OLD DROPS. 

R 


“RTHE ARGUMENTS ARE THE QUANTITIES GIVEN ABOVE IN THAT ORDER. 
RNOTE THAT THE ZEROTH WORD OF THE VECTOR ARGUMENTS ARE TO BE 
-КСІМЕМ» FURTHER THESE VECTORS CONTAIN REQUIRED INPUT AS WELL 
RAS SPACE TO PUT THE DESIRED OUTPUT. FURTHER THE USUAL 
 RPROGRAM COMMON PARAMETERS ARE ASSUMED TO CORRESPOND TO THE 
RARGUMENT TEMPERATURE» 
R 
EXTERNAL FUNCTION (X»DELX»A,»DELA»P»RHO,T,U»M»MDOT»ALPHA »N» 
1 RADsDELR»NDOT»DELG ) 
R 


AP —— —  — — — w — — ———  — 8 








PROGRAM COMMON GAMMA» MU» B» C» CP» L» RHOL» SIGMA 
INTEGER N 


ENTRY TO NUCLE 20 


EXECUTE NUCLEle (A»DELX»MDOT»P>T>RADIN)»NDOT(N)>DELG) 


R WE HAVE NOW ONLY TO COMPUTE THE RADIAL INCREMENTS AND 
К THEIR CONTRIBUTION TO THE CONDENSATE MASS INCREASE» THE 
R DROP TEMPERATURE IS TAKEN AS THE SATURATION CURVE 


R TEMPERATURE CORRESPONDING TO THE GIVEN PRESSURE» 


— = =: 
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К 
TORCP=C/(B-ELOGe (P)) 


“DELR= 7978.846 * SQRTel 1¢379%e6027/(T*MU)) * ( TDROP-T ) i 


DELR=DEL X* (ALPHA/L)*(P/RHOL) *(DELR/U) — 


— THROUGH SI,.FOR Is1ls»1»leGeUN-1) 


EI DELR(1)=DELR an. 
— ИБ > LEM URL. zz DES en ы За”. а M a — t nB HÇ. G -FWm In À————— É—Á 
R____ COMPUTATION OF RESULTANT TOTAL CONDENSATE MASS INCREASE 
R | 
_DEL62=0, ЕЕ 
THROUGH S2.FOR = ОБЕСИ 
52 DELG2-DELG2« RAD(I)*RAD(CI)*NDOT(I)*DELR(I) > 
- DELG=DELG+12+5663706*RHOL*DELG2/MDOT 
R лақа ++ m SEELE SEP. A ro 6 Ar 
FUNCTION RETURN 
R 
INTEGER I 
R 


END OF FUNCTION 





meme PIL’ 
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O CEhIJECI sPUNCH Cor Gi NOZEEQOINGZE COO] 
К AERONAUTICAL ENGINEERING 
К NOZZCE CONDENSATION AREA 
R 
R NOZZEE 
R 
RRESEARCH-- Je GRIFFIN PROGRAMMING== Leo HARDING ОЕТ 2 
R 
RORTGIN OF THE COORDINATE SYSTEM FOR THE МОХА ВЕСЕЛЫЕ ЕЕ Е 2 
ROF ТИЕ ТНКОАТ» І.Е.» А(О) # А5ТАК» THE МАКТАП О БЕЛЕК ТМ Esas n 


RPOSITIVELY IN THE DIRECTION GF FLOW AND LS N ГУЕ КШ 
RFRCM THE THROAT. THE NOZZUE PARAMETERS ARE 


R ASTAR = THRCAT AREA 

К XMIN = X-COORDINATE Cr INTAKE 

R ХМАХ = X-COCRDINATE OF EXIT 

R LIANG = INTAKE HALE ANGEE TNEDEGREES 

R QOJTANG = FXIT STOE THA F ANGLE INED OFRER 

THE EFEANTRY THROAT READS AND PRINTS THESF PARAMETERSS 
R 


EXTERNAL ҒИМСТІОМ (А1»А2»А3З) 
R 
ENTRY TO NOZZLEe 


WHENEVER Al eEo STHROATS 


R 
R RETURN THROAT AREA AND INITIAL UZE SIE ONE GGeooan 
READ DATA 

PRINT COMMENT 50 NOZZLE PARAMETERSS 


PRINT RESULTS ASTARsINANGsOUTANG 
PRINT RESULTS XMINs»XMAX 
A=INANG/RADIAN 
INTAM=SIN.(A)/COS. 14) 
A-OUTANG/RADIAN 
OUTTAN=SIN+(A)/COS.+(4) 


A3=ASTAR 
FUNCTION RETURN 
R 
OR WHENEVER Al eEe SAREAS 
R 
R COMPUTE THE AREA AT A2 AND RETURN IN A323 


WHENEVER (XMINeLeA2) «AND. (А2.(.0.) 
A=INTAN#(.ARS. A2) 
OR WHENEVER (А2.С6.0»ө) «AND. (A2eL eXMAX ) 
A-OUTTAN*A?2 
OTHERWISE 
2-0. 
WHENEVER A2eF e0+«sA3=ASTAR 
FUNCTION RETURN 
ENO OFFC ONDI TIONAL 
A3=ASTAR+20*A 
FUNETTON "RE TURN 


OR WHENEVER Al eFe SINVERSS 


COMPUTE THAT POINT ON THE SANS ORTA ANA Т 
HAS AREA 'A3'. THIS POINT DEPFNDS ON ЕНЕ МАСЕ NUMBER 
WHICH MUST BE IN A2 ON ENTRYs SEE THE PROGRAM БОК 

TRE PARTICULARS se IF A2(1 І5 ЖО RETURN eae 
SOLUTION 1S IN A2s IF Ti 91s ZERO NAO ЕЕ 


AI aX 





po 


cdm STORED IN A2 BUT DOES NOT LIE INSIDE THE NOZZLE. 
MEE AS ASTAR) a 
WHENEVER Ne lene Ol es We Ue 


ас Ет ac NN » i 
MACH=A2 
_ WHENEVER MACHeLole P D! ке. 





Х--А/ІМТАМ 

WHENEVER XMINe«GeXsX(1)70. 
—BTHEBWISE UO mc ue 

X=RZOUTTAN 
0 WHENEVER Х.6.ХМАХЭХ(1150.00 

END OF CONDITIONAL 

EX 
A2(1)=X(1) 


— 
=- 
m a a Rn rn n. saru == 2 `. = 4 TS Ee o 


UNCTION RETURN ` 


+ ere eer ee mmm - 


Saye sas scc P — HP — e — —— a ee el eee 


* e » Фатмақеа, жағ 


R 
END OF CONDITIONAL 
R 
INTEGER Al 
ы ONE) C 
VECTOR VALUES RADIANS57.2951795 
R 


END OF FUNCTION 
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BCOMPILE MADsPRINT ORJECTsPUNCH OBJECT VAPOROOLVAPOROO* 
к AERONAUTICAL ENGINEERING ae N 
R NOZZLE CONDENSATION RESEARCH 
R 
R ~ VAPOR PARAMETERS 
R 
RRESEARCH-- Je GRIFFIN PROGRAMMING-- Le HARDING JULY 1963 
R 
RCOMPUTE AS A FUNCTION OF TEMPERATURE THE FUNCTIONS 
RO RHOL = LIQUID DENSITY — 

R L = LATENT HEAT 
R SIGMA = SURFACE TENSION 


RNOTE THAT EACH OF THE QUANTITIES MAY BE COMPUTED SINGLY AND 
RTHAT THEIR VALUES ARE STORED IN LOW CORE. THIS PERMITS THEM _ 
RTO BE HELD CONSTANT WITHOUT INITIALIZING THESE SUBROUTINES 
“RANY TONE» OR ALL» OF THESE PARAMETERS MAY BE HELD CONSTANT BY 
RSIMPLY RETURNING THE VALUE FROM PROGRAM COMMON. 
RIT IS TO BE NOTED THAT SIGMA AND RHOL ARE FUNCTIONS OF THE 
RDROP TEMPERATURE» NOT THE VAPOR TEMPERATURE.» IT IS THE DROP 
КТЕПИЕКАТУНЕ ТНЕМ THAT MUST BE USED” IN CALLING ON II F eu Sh 
“RENTRIES»s WHEREAS THE VAPOR TEMPERATURE IS USED WHEN COMPUTING 
RTHE LATENT HEAT. 
ur a ne d 
MENA FUNCTION (T) 
б mc S 
PROGRAM COMMON GAMMA» MU» B» C» CP, 1» RHOL» SIGMA 
Әсе” — займа e AAN SA 
ENTRY TO FRHOLe 
R COMPUTE THE VALUE OF THE LIQUID DENSITY 
FUNCTION RETURN VALUE 
Кк | E 
ЖЕСЕ TO P D, м. 
Б COMPUTE THE VALUE O 
FUNCTION RETURN VALUE 
"ODE E 
ENTRY TO FSIGMA. 
R COMPUTE THE VALUE OF THE SURFACE TENSION 
FUNCTION RETURN VALUE _ 
R 
END OF FUNCTION 


ae => AU AM AA A DI s - m s ee eee 70 


LATENT HEAT 





THE 


A m ° , = oP OES SEE sie PE AT Al > DS 





APPENDIX C 


INPUT DATA 


Summary of input data and units required. 

Table of input data used to obtain results plotted in Figures 1 
through 13. 

Plots of vapor properties obtained from References 20 through 
24, 


2592 





Summary of Input Data Requirements 


Parameter Definition Units 
PZERO Initial pressure (Р) dyne/ cm? 
RHZERO Initial density (р) gm/ om 
TZERO Initial temperature (T) | OK 
PE Initial approximation to saturation pres- dyne/ eue 

sure (Р) 
GAMMA Ratio of specific heats (y) N. D. 
MU Molecular weight of the vapor (py) gm/gmol 
CP Specific heat (C) dyne-cm/ gm-?K 
L Latent heat (L) dyne-cm/gm 
SIGMA Surface tension (0) dyne/cm 
RHOL Liquid density (oy) gm/ сш : 
4 ) Saturation curve constants b: 
ALPHA Accommodation coefficient (a) N. D. 
DELX Incremental step in X (Ax) cm 


XRANGE Length of interval, starting with conden- cm 
sation point, over which flow conditions 
will be computed 


XPOINT Points which will be printed out, i.e., N. D. 
every fifth, tenth, etc. , point. 


TRANGE Temperature range over which program ЧК 
will search for condensation point 


DELTAT Initial step size in condensation point OK 
search over TRANGE 


PERCNT Value of epsilon for determination — 


condensation raction 
ASTAR Throat area (A*) cm“ 
AMIN X-coordinate of intake cm 
XMAX X-coordinate of exit cm 
INANG Inlet half angle degrees 
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Summary of Input Data Requirements (continued) 


OUTANG Diffuser half angle (0 1 /9) degrees 

NA Avogadro's number (in program) molecules/gmol 
6, 027 х 1023 

R Universal gas constant (in program) . dyne-cm/ gmol^K 
8.314 x 107 

K Boltzmann's constant (in program) dyne-cm/ ai 


1.379 x 10-16 
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Table C-1. Input Data 


N, (1) N,(II) Cu Zn(1) Zn(II) 
6 
P. 8, 45 x 10 8. 45 x 106 1x 10° 1x 10° 1х 10° 
p. . 00952 , 00961 2. 5A x 1074 3,93 x 1074 6. 04 x 1074 
E 298 295 3000 1500 1300 
p 2000 2000 166. 7 166. 7 166. 7 
Y 1.4 1.4 1. 667 1. 667 1. 667 
n 28. 02 28. 02 63. 54 65. 37 65. 37 
C 1.07 x 10! 1.07 х107 3.26x 10° 3.18x10° 3.18 x 10° 
p 9 9 10 10 10 
ү; 2.3 x 10 2.3 x 10 4. 94 x 10 1. 82 x 10 1. 82 x 10 
g 25685, 297 24, 25-. 22T 1300 1000-. 267 1000-. 26T 
Pr, 1. 177-.00476T 1. 177-. 00476T 8.9-. 0007T 7.4-.0008T 7.4-. 0008 
С 882 882 34, 200 13, 800 13, 800 
B 25.7 25.7 25. 46 25.3 25.3 
о 1.0 1.0 1.0 1.0 1.0 
А* . 0645 ‚ 0645 „ 645 ‚ 645 ‚ 645 
<a -5 -5 -5 -5 -5 
min 
25 25 58 21.1 58 
max 
INANG 45 45 45 45 45 


OUTANG 9. 75 Je 20 45 20 





2) 


P (dynes/cm 


10 


10 


10 


mM 
сә 
bo 


10 | 


1 
3 


E === 
o EEE ====== 
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EXE SERRE Au =: "C-C Curve Match 
= ——— a 


EE == =1np = 25.7 - 882/T К 
со йк 


== 
== == == =: E E = киси 
== = = p x == Matched at 40° and 60° 
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Figure C-1. Vapor Pressure of Nitrogen,NBS Circular 564 (24) 
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Figure C-3. 
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Figure C-4. Metal Vapor Properties 
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Figure C-5. 






































